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1 INTRODUCTION 


Because they provide high stiffness and strength and are relatively light 
weight, continuous fiber reinforced laminated composite materials are candi- 
dates for use in large space structures. One of the significant differences 
between these composites and more conventional engineering materials is that 
laminated composites develop a considerable amount of large-scale damage. 
Numerous ply cracks and delaminations may develop in a laminate early in its 
load history. 

Unlike fiber breaks and fiber-matrix debonds, which tend to influence 
only a small volume of material, ply cracks and delaminations can cause sig- 
nificant changes in the bulk properties of a laminate. Numerous researchers 
have studied and modelled the influence of ply cracking [1,2,5] and 
delamination [3,6] on the elastic behavior of laminated composites. But while 
significant increases in damping have been observed as a result of damage 
development [7,8], and modelling efforts have predicted that significant 
changes in the dynamic response of structures can result from changes in mate- 
rial damping [9,10], existing models of damping in composites [11,12] have 
not included the effects of damage on damping. 

This report summarizes the work completed under NASA Grant NAG- 9 -192. 

This research effort had as its objective the development of a damage depen- 
dent constitutive model for predicting the influence of damage on the damping 
properties of laminated composites. In particular, the beam bending behavior 
of cross-ply laminates was studied. The damage mode of interest was ply 
cracking. The research proceeded along two fronts. The primary effort was to 
develop a continuum damage model based on internal state variables to describe 
the viscoelastic behavior of laminates as a function ply cracking. In support 
of the continuum modelling, an experimental program was implemented to deter- 
mine stiffness and damping as a function of damage. A second effort involved 
the development of a micromechanics model of damage dependent damping. 
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2 CONTINUUM DAMAGE MODEL 


2.1 Theoretical Development 

The constitutive model used in this study follows the model developed [4] 
and used by Allen, et al . to predict the damage dependent elastic response of 
laminated composites containing ply cracks [5] and delaminations [6]. For a 
given ply, the uniaxial constitutive equation has the form 

• ( 2 * 1 ) 

where E is the modulus of the material, and a! is an internal state variable 
which reflects the effect of ply cracking on ply behavior. In general, there 
would be one internal state variable (ISV) for each different damage mode. 
Multiaxial behavior requires a tensor quantity to characterize the effect of 
damage . 

In the earliest version of the model [13,14,15], it was assumed that a 
ply could be modelled as a Voigt material, which has the mechanical analog 
shown in Fig. 1. For such a material, the constitutive equation has the form 

°xx“^ e xx +T l G xx’ (2*2) 

where the term r\e xx represents the viscous part of the material response. 

This equation can be made to fit the form Eqn. (1) by taking a l = -je xx . In 
fact, two IS V's, one representing stiffness changes due to ply cracking and 
the other representing damping changes due to ply cracking, were used in the 
analysis. Strictly speaking continuous fiber reinforced materials do not 
behave like Voigt materials, within a limited range of frequencies and ampli- 
tudes appropriate Voigt material parameters can be chosen to provide a reason- 
able representation of ply response. 
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While the Voigt model was useful in the early stages of the development 
of the damage dependent constitutive model, it was somewhat restrictive. For 
example, Voigt type ply response led to a linear frequency dependence on lami- 
nate damping [13,14,15]. This frequency dependence was not observed in subse- 
quent experimental results. A more general formulation has since been 
developed which uses a complex modulus representation of viscoelastic ply 
properties [16]. It is this second formulation which will be discussed here. 
In this formulation, the damage dependent, elastic constitutive relation has 
the form 


+ - C 2 - 3 ) 

where now the term la, captures the damage dependence. For a linear viscoe- 
lastic material under sinusoidal loading, the constitutive relation takes the 
form 




* ** 
a, 


(2.4) 


where 


£’-£' + i£ 




a 
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Using the viscoelastic constitutive relation above, it can be shown that 
the energy dissipated by a ply per cycle of loading is given by 

Aw - Jt£'e* + nI'a lA £ A , (2.5) 

and the stored energy during the peak displacement is given by 

U + re A a lA }/2 . (2.6) 

It then follows that the loss factor is 


£'e* + /' a lB 6, 


With the additional assumptions that 


I m a 


IA 




/'a 


1 A 


A£'e„ 


(2.7) 


( 2 . 8 ) 


we arrive at the expression for the loss factor 


T] - 


F- + AE- 
£'+A£‘ * 


(2.9) 


The complex modulus of a laminate may be determined from the ply level 
constitutive relations using classical lamination theory [17] to be 


£* “ 



Zl-lWi 


( 2 . 10 ) 


where E * is the complex modulus of the laminate, h is the thickness of the 
laminate, N is the number of plies in the laminate, E * is the damage dependent 
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complex modulus of ply i, and the z ( are the coordinates of the ply inter- 
faces. For a symmetric cross -ply laminate, a symmetric laminate consisting of 
0° and 90° plies, the expression above reduces to 

E'-aE' L +bE' T (2.11) 

if the laminate is undamaged. In Eqn. (2.11), a and b depend on laminate 
geometry, and E\ and E* T are the complex moduli of the 0° and 90° plies, respec- 
tively. The flexural damping is given by 


tan - r| - 


aE~ L +bE~ 7 


( 2 . 12 ) 


aE\+bE'j 

In the present study, it was assumed that ply cracking would affect E\ 

but not E[ . It follows that for a damaged laminate, the flexural damping is 
given by 


aE' £ +b£' T + CAE' T 


(2.13) 


' a£' £ +6£' T + CA£' T 

where CA£" T and CAf ' r are functions of the total crack surface area and loca- 
tion of the cracks. Both A £' 7 and A E' T can be determined from experimental 
data from a single laminate, and then can be used to predict behavior in other 
laminates via Eqn. (2.13). 


2.2 Experimental Program 


In order to evaluate the analytical model, damping was measured as a 
function of damage in a variety of graphite/epoxy cross-ply laminates [14,15]. 
Damage was introduced into straight sided coupon type specimens under uniaxial 
tensile loading in an MTS servohydraulic testing machine. Edge replication 
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was used to monitor the development of 90° ply cracks during tensile loading. 
At various stages of damage development, the tensile test was interrupted, so 
that the specimen could be removed and its flexural damping measured. 

The experimental setup used in performing the damping measurements is 
shown schematically in Fig. 2. The specimen was cantilevered from a support 
housed in a vacuum chamber. All damping data was obtained in vacuum. A 
motorized wheel with a contacting rod was used to deflect the specimen and 
excite its first mode of free vibration. Care was used to insure that the 
initial deflection was small, so that no new damage would be introduced during 
the damping test. A strain gage mounted on the specimen near the cantilevered 
support was used to monitor the surface strain in the specimen during free 
vibration. A personal computer with an A/D board was used to start and stop 
the motor, and to record strain data. 

After capturing the strain data, the peaks in the strain versus time 
curve were obtained. A representative plot of log strain amplitude versus 
cycle number is shown in Fig. 3. The decay of strain amplitude is seen to be 
nearly linear on a semi-log plot. The slope of the best fit straight line to 
the log strain amplitude data was taken to be the logarithmic decrement, which 
is directly related to damping. 
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Figure 2. Schematic diagram of the damping measurement sys- 






2.3 Comparison of Theory and Experiment 


From measured values of the stiffness and damping of laminates containing 
all 0° and all 90° plies, the parameters E' Lt E' T , E" u and E~ T can be deter- 
mined. These parameters can in turn be used to predict the damping of other 
laminates consisting of combinations of 0° and 90° plies. A comparison of the 
predicted and observed damping in several cross -ply laminates is presented in 
Fig. 4. Agreement between theory and experiment is quite reasonable. The 90° 
plies exhibit matrix dominated behavior, and hence are the primary contrib- 
utors to damping in laminates. In addition, the effect of a 90° ply increases 
with distance from the midplane. 

In order to predict the influence of damage on damping, it was necessary 
to determine the parameters A E' r and A £' T as a function of crack surface area. 
Since the elastic properties of the graphite/epoxy laminates are dominated by 
the 0° plies, A E' T term was assumed to have a negligible effect on the denomi- 
nator of the right hand side of Eqn. (2.13). Further, damage dependent damp- 
ing data from [ 0/90/0] s specimens was used to determine 

A£~ T (s) = 57.97x 1 0~ 6 [0.0002(s) - 1 .7108 x 10’ 5 (s 2 )] , 

where s is the crack surface area in the 90° ply. Using this expression in 
Eqn, (2.13), damping was predicted as a function of damage in other cross-ply 
laminates . 

Figures 5, 6, and 7 present a comparison between predicted and observed 
damping values in [ 90/0/90] Sf [0/90] 2s , and [90/0/90 2 ] s laminates at two dif- 
ferent damage states. Generally, agreement between theory and experiment is 
quite good. The effect of damage on the behavior of the [ 0/90 ] 2s is smaller 
than the effect of damage on the behavior of the other laminates because the 
90° plies are located closer to the midplane. Thus, any influence of the 90° 
plies on the flexural damping is reduced. This stacking sequence dependence 
is also reflected in a smaller damping value for the [ 0/90 ] 2s laminate. 
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3 MICROMECHANICS MODEL 


3.1 Theoretical Development 

In order to better understand the details of damage- induced damping, a 
micromechanics model of damaged crossply laminates was developed [19]. The 
model is based on shear lag theory, and explicitly includes 90° ply cracks. 

In fact, two versions of the model were developed. In the first, each ply is 
treated as an Euler-Bernoulli beam, while in the second, each ply is treated 
as a Timoshenko beam. The Timoshenko beam version will be described below. 

For the Timoshenko model, the axial displacement in layer i is given by 

uXx.zJ- u 0 Xx) + zk,(x) . (3.1) 

where 

• (3-2) 

In Eqn. (3.2) w is the deflection of the ply, which is assumed to be the same 
as the deflection of the laminate. The Euler-Bernoulli model does not include 
the ply shear term, which is represented by [3. 

Ply equilibrium equations are obtained by performing force and moment 
summations on a differential element taken from the ply. Such an element is 
shown in Fig. 8. The interlaminar shear stresses shown in this figure are 
determined using the shear lag assumption, i.e., they are assumed to be pro- 
portional to the difference in ply displacements across the interface. 

In fact, the shear strain in the interfacial region is given by the differ- 
ence in the axial displacements of the two plies divided by the thickness of 
the shear transfer region. Stresses are obtained by multiplying by an 
appropriate shear modulus for the shear transfer region. Following this pro- 
cedure, the shear lag assumption yields 
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(3.3) 


G( a, a,. i ^ 

**4.4-1 “ ^ I U 0i 2 U °t- 1 2 ^ l ' l J 

_ of a a,, ^ 

u oi 2 u omi “ 2 *** 1 J ’ 

Axial equilibrium for the ply has the form 


G f a<-i a,., , * - 

— * 4M j U„ -0 

equilibrium in the z direction yields 


(3.4) 


°«,' = 0 . 


(3.5) 


and moment equilibrium has the form 


6 G ( , a,., a,.. 


1 2G , 




?k t 



kS 


\2G X . V 

A - , 

J-l 


(3.6) 


Exterior plies have slightly different equilibrium equations owing to the fact 
that on one side of such a ply there is no neighboring ply, and thus there are 
no interlaminar shear stresses. 

The equilibrium equations for the plies within a laminate can be 
assembled into a system of coupled, second order differential equations in the 
form 


£{y}-{y~> - {c> 


(3.7) 
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where the column vector { y } has as its components the midply displacements u 0 i 
and the ply rotations k im These governing equations are applied to a region 
of the laminate which corresponds to one -half of the crack spacing. Such a 
region is shown in Fig. 9. At x-0, the displacement in the uncracked plies is 
zero, while in the cracked plies, the crack surfaces are stress -free. A 
symmetry condition at x-1 requires that ply displacements be equal at the ply 
interfaces. Applied loads then determine the midplane displacement and rota- 
tion at this end. 

The system of linear ordinary differential equations was solved using 
techniques commonly used for analyzing dynamic systems [20]. This approach 
requires the determination of eigenvalues and eigenvectors of the coefficient 
matrix [B] . The solution takes the form of a linear combination of "mode 
shapes" : 


{y}“X“ ln/1 < e° ,X <V ( } . (3-8) 

1 

where the a,- are the eigenvalues and the {ip*} are the eigenvectors. Boundary 
conditions determine the coefficients A i in the solution. 

The viscoelastic problem was also solved in order to determine the effect 
of cracking on damping. This was accomplished using complex moduli in the 
analysis described above, and obtaining the corresponding solution. Note that 
only material viscoelastic behavior was included, and that other dissipation 
mechanisms, such as friction, were ignored. 
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3.2 Comparison with Experiment 

Experimental data for bending stiffness as a function of crack density 
was obtained for [0/90 3 ] s and [90 2 /0 2 ] B glass/epoxy laminates. A glass/epoxy 
system was chosen because it is a less fiber dominated system than graphite/e- 
poxy, and thus should exhibit larger property changes due to damage develop- 
ment. Specimens were loaded in quasi-static tension in order to introduce 
damage, which was monitored using edge replication. Four-point bending tests 
were used to determine the bending stiffness at a variety of damage states. 
This data was compared to model predictions. 

Figures 10 and 11 present the theoretical and experimental bending stiff- 
ness data for the [0/90 3 ] s and [90 2 /0 2 ] s laminates, respectively. Note that in 
Fig. 10, there is little bending stiffness reduction observed. This is 
because the 0° plies are located at the exterior of the laminate. In light of 
this small stiffness reduction, the predicted responses are in reasonable 
agreement with experiment. The Timoshenko model provides slightly better pre- 
dictions than the Euler -Bernouli model. The [90 2 /0 2 ] s laminate (Fig. 11) 
exhibits larger stiffness changes. The Timoshenko model provides a reasonable 
estimate of laminate response up to crack densities of about 13 cracks/in. , at 
which point delaminations begin to develop. 

Fig. 12 shows the predictions of loss factor as a function of crack den- 
sity in a [90 2 /0 2 ] s laminate. At low crack densities, the damping increases 
with cracking. But at larger crack densities, damping decreases with 
cracking. It is believed that this behavior is a result of excessive load 
transfer from the relatively viscoelastic 90° plies to the relatively elastic 
0° plies that occurs at high crack densities. At lower crack densities, 
stress concentrations in the shear transfer region are able to generate an 
increase in damping. This decrease in damping is not observed in the labora- 
tory. At present, it is believed that other dissipation mechanisms , such as 
friction between rubbing crack faces, account for the discrepancy. With 
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Normalized Bending Stiffness 



Figure 10. Damage dependent bending stiffness in [0/90 3 ] s 
glass/epoxy laminates. 
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additional effort, it should be possible to use the micromechanics model to 
estimate frictional dissipation. 
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Figure 12 . Predictions of loss factor as a function of damage 
using the Timoshenko model. 
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4 CONCLUSIONS 


A constitutive model has been developed for predicting damping as a func- 
tion of damage in continuous fiber reinforced laminated composites. The dam- 
age model is a continuum formulation, and uses internal state variables to 
quantify damage and its subsequent effect on material response. The model is 
sensitive to the stacking sequence of the laminate. Given appropriate base- 
line data from unidirectional material, and damping as a function of damage in 
one crossply laminate, damping can be predicted as a function of damage in 
other crossply laminates. Agreement between theory and experiment has been 
quite good. 

A micromechanics model has also been developed for examining the influ- 
ence of damage on damping. This model explicitly includes crack surfaces. 

The model provides reasonable predictions of bending stiffness as a function 
of damage. Damping predictions are not in agreement with experiment. This is 
thought to be a result of dissipation mechanisms such as friction, which are 
not presently included in the analysis. 
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ABSTRACT 


. An Experimental Investigation of Damage-Dependent Material Damping 
of Laminated Composites. (May 198S) 

Scott Andrew Smith, B.S., Texas A&M University 

Co-Chairs of Advisory Committee: Dr. Alton L. Highsmith 

Dr. Charles E. Harris 

An experimental program was developed in which the material damping 
of laminated graphite /epoxy specimens was measured as a function of matrix 
cracking. The damping test apparatus was designed to measure the first mode 
free vibration response.of a cantilevered beam. Matrix cracks were introduced 
into cross-ply laminates, and damping, calculated by the logarithmic 
decrement method, was measured at several levels of damage. 

The 90 degree plies were found to have a greater damping capacity than 
the 0 degree plies, and placement of the 90 degree plies away from the mid- 
plane of the laminate thickness resulted in a higher damping measurement. 
Matrix cracks in the 90 degree plies were seen to significantly increase the 
damping of the cross-ply laminates. In addition, the matrix cracks which 
formed in 90 degree plies near the specimen surface, caused a greater increase 
in damping than matrix cracks near the mid-plane. Finally, a damage- 
dependent material damping model was applied, which closely predicted the 
measured damping values. 


31 



ABSTRACT 


Matrix Cracking and Bending Stiffness Reduction 
in Composite Laminates. (December 1988) 

James Alan Frailev, B.S., Texas A&M University 
Chair of Advisory Committee: Dr. Alton L. Highsmith 

Shear-lag theory was used in the development of two mathematical models 
to predict the behavior of cracked, cross-ply, laminated composites. One model used 
Euler-Bernoulli beam theory to model the deformation of each ply while the other 
used the Timoshenko beam theory. Solutions to the governing system of differential 
equations provide displacements and stresses throughout a representative volume 
of the laminate. The models were used to predict the response of laminates under 
two loading conditions: axial, pure bending. 

Analytical results were compared to experimental data for axial and bend- 
ing stiffness reductions with matrix crack density. Glass/epoxy [O/ 9 O 3 ] j and 
teOj/Ojl* laminates were selected for the analytical— experimental comparison. 
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ABSTRACT 


An Investigation of Damping Behavior of 
Fiber Reinforced Composite Materials. (May 1989) 

Shankar Kalvanasundaram, B.Tech., Madras Institute of Technology 
M.S., Oklahoma State University 

Co-Chairmen of Advisory Committee: Dr. David H. Allen 

Dr. Walter E. Haisler 

Advanced fibrous composite materials are known to exhibit viscoelastic be- 
havior, which in turn contributes to the damping of these materials. This study 
has investigated the effect of shear deformations, rotary inertia and matrix crack 
damage on the damping characteristics of composite materials. For a composite 
beam the influence of shear deformation and rotary inertia depend on the ratios of 
E' /G' , T an /Tan§ E , the length to thickness ratio and mode number. For highly 
anisotropic beams the loss factors in shear are substantially higher than longitudinal 
loss factor. This leads to large differences between the damping predictions obtained 
from Euler-Bernoulli and Timoshenko beam theories. The results of this study show 
that for highly anisotropic slender beams in which the ratio of Tan$ q /T an® £ and 
E’ /G‘ are high, the influence of shear damping is important even in low vibrational 
modes. A constitutive model has been developed for predicting the effect of mi- 
crostructural damage on the dynamic response of composite materials. This model 
is valid for a fixed damage state and small amplitude vibrations. The theoretical 
formulation has been verified for crossply graphite-epoxy laminates with transverse 
matrix crack damage. For crossply laminates the loss factor in shear is of the same 
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order as the longitudinal loss factor and Euler- Bernoulli theory has been found to 
be adequate for obtaining dynamic properties. Damage dependent material con- 
stants obtained from [0/90/0], laminates have been used to predict the increase 
in damping in [90/0/90],, [0/90j 2 , and [90/0/90 2 ], laminates for different matrix 
crack damage states. Damping has been found to be more sensitive to matrix crack 
damage than the stifness loss. Thus, damping measurements hold some promise for 
further studying the damage development in composite structural components. 


34 


8 TECHNICAL REPORTS 


35 


8.1 


Dynamic Response of a Viscoelastic Timoshenko Beam 



DYNAMIC RESPONSE OF A VISCOELASTIC TIMOSHENKO BEAM 

by 

S. Kalyanasundaram* 

D.H. Allen** 
and 

R.A. Schapery*** 

Texas A&M University 
College Station, TX 77843 


Abstract 


Nomenclature 


The analytical determination of the* steady 
state and transient flexural vibrations of beams 
Is often based on the Euler-Bernoulll equation. 
In this theory what are normally secondary 
effects, such as shear deformations and rotary 
inertia are not included; but these effects may 
significantly influence the vibrational response 
under certain conditions. 

Advanced fibrous composite materials are 
known to exhibit viscoelastic behavior, which in 
turn contributes to the damping of these 

materials. The analysis presented in this study 
deals with the vibratory response of viscoelastic 
Timoshenko beams under the assumption of small 
material loss tangents. The appropriate method 
of analysis employed here may be applied to more 
complex structures. 

This study compares the damping ratios 
obtained from the Timoshenko and Euler-Bernoul 11 
theories for a given viscoelastic material 
system. From this study the effect of shear 
deformation and rotary inertia on damping ratios 
can be identified. 
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Introduction 


The study of vibration damping properties of 
advanced composites is important for several 
reasons. An automatically controlled flexible 
space structure may act either docile or 
unmanageable depending on the level of damping in 
higher vibration modes*. It is also known that 
damping is sensitive to microstructural detail so 
that damping measurements can be used to study 
damage development in composite members . 

The analytical determination of the dynamic 
properties or beams generally employs the Euler- 
Bernoulli equation. For highly anisotropic 
beams, neglecting shear deformations and rotary 


1 



inertia can result in relatively large errors In 
the predicted dynamic properties. For long 
slender beams composed of undirectional composite 
laminae in which the ratio of longitudinal 
modulus (£j|) to the shear modulus (G 12 ) Is high, 
transverse shear deformations can become 
significant even in lower vibrational modes. 

In experiments^ conducted on unidirectional 
E-glass composite beams, the measured values of 
longitudinal loss modulus are typically 50 to 80 
percent higher than the prediction obtained from 
Euler-Bernoull i beam theory. Adams, et al. 4 " 5 
have also indicated that the measured damping 

values may depend on shear effects. Schultz and 

7 8 

Tsai have used the damping ratios measured for 
three fiber orientation angles to predict the 
damping ratio for any angle. The damping ratios 
exhibited qualitative but not quantitative 
agreement with the predicted properties. 

There is a need to develop an improved model 
for the analysis of large space structures which 
is capable of describing global structural 
response including accurate representation of 
damping characteristics of structural 
components. The global dynamic response of space 
structures can be modelled with equivalent 

Q 

continuum beam and plate models. Sun, et al." 
have presented procedures for constructing 
several equivalent continuum models. It was 
found that tne Timoshenko beam and Mindlin plate 
theories were suitable for equivalent continuum 
modelling. Material damping 1 ^ of structural 
members can be included in such formulations. 
For small material loss tangents the dynamic 
steady-state and transient response of a 
viscoelastic Timoshenko beam and other structures 
can be predicted by simple techniques described 
by Schapery 11 . 

The Timoshenko beam theory 1 ^, which includes 
the effects of shear deformation and rotary 
inertia, has been used to predict flexural 
resonance frequencies obtained for anisotropic 
specimens such as wood 1 ^ and fiber reinforced 

i * 

composites Although the Timoshenko theory was 
originally oeveloDed for isotropic materials, it 
may be applied to anisotropic beams, providing 
the appropriate longitudinal and shear moduli and 
associated loss factors are employed in the 


model. When the wave lengths of the deformation 
response are relatively small, a higher order 
shear deformable theory may be necessary In order 
to accurately predict dynamic properties of 
composite beams 15,16 . 

The goal of the present investigation Is to 
compare the f ree-vibration damping obtained from 
Timoshenko and Euler-Bernou 1 1 i beam theories for 
a simply supported beam. From this study the 
effect of shear deformation and rotary inertia on 
damping can be identified. The results of the 
study can be used to interpret the difference 
between experimental and predicted damping values 
for beams and to guide a more general study of 
passive damping in large space structures. 


Analysis 

The governing differential equation of the 
free vibrations of an elastic EOler-Bernoul 1i 
beam is given by 


4 2 

EI — + oA -Mr = 0 

SX St 


( 1 ) 


The solution for free vibrations can be written 
in the form 

w(x,t) = W(x) exp(iut) (2) 


where W(x) is the displacement amplitude at any 
point in the beam. For sufficiently small 
damping, free vibrations of a viscoelastic 
structure decay very slowly, which permits the 
use of an approximation wherein stress-strain 
equations for steady-state harmonic vibration are 
employed. Thus, by the correspondence principle 
of linear viscoelasticity, the elastic Young's 
modulus E is replaced by the corresponding 
complex modulus E* , where 


E * (1-i Tans.) 


( 3 ) 
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which. In general, depends on frequency. The 
boundary conditions for the simply supported 
Eu ler-Bernoull i beam are 


w 


* 2 W 

“I? 


0 


at x*0, L 


(«) 


The deflection of the nth mode of vibration is 
given by 


w n (x.t) = A 


n 


sin 


nnx 

L 


Jut 
e n 


(n-1.2....) (5) 


Tan*, 


'(I - 


2 

n n 

/lJ n } w n 


( 9 ) 


which shows that w'/u' is small when the 
n n 

material loss tangent tan*^ Is small. 

The coupled equations for the total 

deflection w and bending slope for the 

Timoshenko Deam^ are as follows: 


.2 


El -2-5 + k (77 - *) AG - Ip ■ 2 -y * 0 (10) 

bx^ oX at 4 


The frequency equation is given by 
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( 6 ) 


Since E* is complex, the frequency values are 
also complex and can be written as 


w + 1 u> 
n n 


( 7 ) 


Eliminating * from equations (10) and (11), the 
following equation in w is obtained: 


n + i!- _ | (1+E/kG) 2jfo 
p * sx* 3 c M ax*at* 


♦kit !- 0 

O C 


( 12 ) 


t i 

Observe that when u> n > 0, the displacement, 
equation (5), is predicted to decay 
exponentially. Separating the real and imaginary 
parts of equation (6) results in 


E = 


dAL h '2 n w n 

4 4. w n ( 1 " 


nrrl 


if k'/O 2 «i. 


Tne boundary conditions for the simply supported 
Timoshenko beam are given by 


w = - 0 at x-0,1 

oX 


(13) 


The deflection of the nth mode of vibration is 
given by 


dAL 

4 4 T 
n n I 
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( 8 ) 


w n (x,t) 


A sin 
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mrx 
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Jut 
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(14) 


As in the Euler-Bernoul 1 i theory for small 
Also, damping, tne complex frequency equation is 

obtained by using equation (14) in (12) and by 
replacing E and G by tne complex moduli E and 
G*, where 
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E - £*(1 -r i Tan* E ) 


and 


(15) or 
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Therefore, 


t I n n , , n r I / 1 , t \ i Z 

■" T ^ >l ‘» 


fcl 

kG*A 


(17) 


whicn is identical to the frequency equation for 

» » 

an elastic beam except E and G appear in place 
of elastic constants. Similarly, by setting the 
imaginary part equal to zero, the following 
equation is obtained: 


Now consider again equation (7) and the small 
loss tangent approximation, such that 


•'2 

2 2 

•• 1, Tan 6^ << 1, Tan « 1 
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u n 
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* -^(1 - 41 4 ) (19) 
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Using (19) in (17) and separating the real and 
imaginary parts results in two equations. 
Setting the real part equal to zero gives the 
following equation: 
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I 

Observe that u can be obtained from eouation 
n 

(21) for each n. This value can then be 
substituted into equation (22) to obtain 

i i t 

u . u can be easily estimated 

n n 

2 2 2 

when n « I/L A is a small quantity. Using this 
approximation, it can be shown ** that the last 

term in equation (21) is small compared 

2 2 2 

to . n i I/LA. Neglecting the last term, the 

i 

frequency « n is given by 

22 1 n 2 2 i c' 

-S- l 1 - 7^4-1 I 1 * ^-rlJ < 23 ) 

\r L L kG 



Tanc for the viscoelastic Timoshenko beam is 
defined as. 


4 



2 u 

Tan* - — (24) 

“n 

where J and are found from equation (22). 
n n 

This loss tangent provides a measure of the decay 

per cycle of free vibration, which can be shown 

by substituting equation (7) into (5) . and 
• « 

replacing u> n in favor of tan* by means of 

equation (24). The ratio of displacement at 

time t + 2«/u' to that at time t Is found to 
n 

be exp(-ntanc) . Note that tan* * tan* £ for the 
Euler-Bernoulli theory, equation (9). 


Discussion of Results 

The influence of shear deformation on the 
damping will depend on the ratios of 

• i 

E /G , Tanc G /Tan*£, the length to thickness 
ratio, and mode number. For unidirectional 

i i 

composites, the ratios E /G and Tan* £ /Tan* G 
depend on the volume fraction of fibers and. 

matrix and the moduli of fibers and matrix. By 

considering different material systems and 
different orientation angles it is possible to 

have wide variations in the ratios 

i i 

E/G and Tan*^/Tan*^. In the continuum model of 
large space structures, the effective shear and 
longitudinal moduli can be varied independently 
of each other by changing the conf iguration of 
the individual cells. All the results presented 
in this study will be given in a dimensionless 
form. 

The material properties and dimensions of 
the beam are given in Table 1. For the purpose 
of illustration the loss factors and storage 
modulus are assumed to be independent of the 
frequency. 


Table 1 

Typical Properties of A Composite Beam 

i 

f’ « 20xl0 6 psi 10-40 

G 

L * 50 In ^ - 20-100 
n 

Tan* r 

Tan *G “ °- 05 Tan*^ * “O' 1 

p e 0.065 lb/in^ width b * 1.0 in 

shape factor k = 5/6 

i < 

The influence of E /G, Tan* G /Tan* £ , L/h and mode 
numbers on the ratio of damping predicted by the 
Timoshenko and Euler-Bernoul 1 i theories are shown 
in Figs. 1 through 8. 

Figure 1 is a plot of Tanc/Tanc £ for 
Tanig/Tanc^ varying from 100 to 1. The first 

mode of vibration (n=l) is considered. The 

length to tnickness ratio of the beam is 100. It 
is found that for large values 
of Tan* G /Tanc^ (100) the damping values predicted 
by Timoshenko theory are 9% to 33% higher than 
the prediction by Euler-Bernoul 1 i theory 

i i 

depending on the ratio of E /G. 

i i 

Increasing E /G increases the shear effects. 

The third mode of vibration (n=3) is 
considered in Fig. 2. The effect of shear 
damping is more pronounced than in the first 

mode. Figures 3 and 4 illustrate similar trends 
for a beam of smaller length to thickness ratio 
(50). The damping Tan* predicted by the 

Timoshenko theory is substantially higher than 
Tan; E ( 38 % to 1300% for Tan* G /Tanc £ = 100). 
It is interesting to note that when the 
difference between Tanc^ and Tan^ becomes small 
the difference in predicted damping values 
obtained from the Timoshenko and Eul er-Bernoul 1 i 
theories is very small. 


In Figs. 5 and 6 damping is predicted for 
different modes while keeping the ratio of 



» I 

E /G constant. For higher modes the shear 
effects become important and this effect Is 
enhanced as the ratio of Tan+g/Tan+g Is 
increased. 

Figures 7 and 8 illustrate the effect of the 
variation of length to thickness ratio on 
damping. For small L/h ratios the damping 
predicted by Timoshenko theory is higher than the 
prediction by Euler-Bernoull 1 beam theory. 

Conclusions 

This study has investigated the effects of 
shear deformations and rotary inertia on the 
damping characteristics of beams. For typical 
longitudinally reinforced composite beams, the 
loss factors in shear (Tancg) . are substantially 
higher than Tan$£. This leads to large 
differences between the predictions obtained from 
Euler-Bernou 1 1 i and Timoshenko beam theory. The 
results of this study indicate that for long 
slender oeams in which the ratio of 
Tan^/Tan^ and E’/G’ is large, the influence of 
shear damping is important even in low 
vibrational modes. The methods and results from 
this investigation can be used to study the role 
of snear damping in higher modes of vibration of 
large space structures. 
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Fig. 1 Effect of sheer deformation on damping es 
a function of Tanc p /Tanc c for various 
values of E*/G ‘ (L7h = 100, n=l) 



Fig. 2 Effect of shear deformation on damping as 
a function of Tana./Tanc- for various 
values of E'/G’ (L/h = 100, n=3) 
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ABSTRACT 

The equations of motion are developed herein 
for laminated composite beams with load-induced 
uatrix cracking. The damage Is accounted for by 
utilizing Internal state variables. The net 
resu.lt of -these variables on the field equations 
Is the introduction of both enhanced damping, and 
degraded stiffness. Both quantities are history 
dependent and spatially variable, thus resulting 
in nonlinear equations of motion. It Is 
explained briefly how these equations may be 
cutsl-llnearized for laminated polymeric 
composites under certain types of structural 
loading. The coupled heat conduction equation Is 
developed, and It is shown that an enhanced Zener 
damping effect Is produced by the Introduction of 
relcrostructural damage. The resulting equations 
ere utilized to demonstrate how damage dependent 
material properties may be obtained from dynamic 
experiments. finally, experimental results are 
compared to model predictions for several 
composite layups. 


E o . - viscoelastic material parameters 

E - axial modulus of elasticity 

- axial modulus of elasticity of 1th 
lamina 

Ej - reference axial modulus of 

x elasticity 

f - transverse load applied to free 

end of cantilever beam 

- transverse load amplitude 

h - beam depth 

I* - modulus weighted moment of inertia 

W about y axis 

k - material constant In Internal 

state variable growth law 

k^j - coefficient of heat conduction 


NOHEKCLATURE 


b 

C 

c i. 


- beam cross-sectional area 

- cross-sectional area of 1th lamina 

- modulus weighted beam cross- 
sectional area 

- coefficients In cantilever beam 
displacement field 

- beam width 

- viscoelastic modulus 

- specific heat at constant volume 

- coefficients In cantilever beam 
displacement field 


*D1j 


-loss In heat conductl on 
coefficient due to damage 

- beam length 

- beam moment resultant about y axis 


y 

N 

n 

n 2 

P 

P T 


- beam moment resultant about y axis 
due to temperature change 

- number of cycles of applied load 

- number of plies In laminate 

- material constant In Internal 

state variable growth law 

- beam axial resultant 

- beam axial resultant due to 

temperature change 
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- beam axial load per unit length 

- bean transverse load per unit 
length 

- heat flux vector 

- heat source 

- tern used to obtain laminate 
damping coefficient 

- surface area of matrix cracks in 
local volume element 

- temperature 

- reference temperature 

- time 

- axial displacement 

- crack opening displacement 

- axial displacement of neutral 
surface 

- Internal energy dissipation per 
unit volume due to matrix cracking 

- total Internal energy dissipation 
due to matrix cracking 

• beam shear resultant In 2 
direction 

- transverse displacement of neutral 
surface 

- beam axial coordinate direction- 

- beam transverse coordinate 
direction 

- 2 component of modulus weighted 
centroidal axis 

- 2 component of distance from 
centroid of 1th ply to modulus 
weighted centroid 

- coefficient of thermal expansion 

- Internal state variable (ISV) due 
to damage 

- axial component of matrix cracking 
damage ISV 

- axial component of matrix cracking 
damage ISV In 1th lamina 

- initial displacement applied to 
free end of cantilever beam 

- axial strain component 


- average axial strain In 1th lamina 

- damping coefficient 


- axial damping coefficient In 1th 
lamina 

- bending damping coefficient In 1th 
lamina 

- axial stiffness loss coefficient 
In 1th lamina 

• bending stiffness loss coefficient 
in 1th lamina 

- modal frequency of mth resonant 
mode 

• mass density 

- stress amplitude of cyclic loading 

- axial stress component 

- shear stress component 

- temperature minus reference 
temperature 

- beam rotation about y axis 


INTRODUCTION 

Mlcrostructural damage In laminated fibrous 
composites has been studied by several 
researchers In the past decade. It has been 
observed that this damage can result In 
significant stiffness loss 11-3). Furthermore, 
experimental results indicate that material 
damping may increase by an order of magnitude or 
more as a result of damage |<,5]. Recent 
research has Indicated that this damage* Induced 
stiffness loss may have a substantial effect on 
the dynamic response of both truss structures [6] 
and beams In bending [7). However, to these 
authors 1 knowledge, no concerted effort has been 
previously reported to model the effect of damage 
induced damping on dynamic structural response. 
While there have been models 18,9] developed for 
damping in laminated composites, these models do 
not include the effects of microstructural 
damage. The primary objective of the research 
described herein Is to develop models that 
explicitly account for the effect of 

mlcrostructural damage on the stiffness and 
damping In the structural dynamics equations of 
motion. 

In this paper a model Is developed for 

predicting the structural dynamics response of 
composite beams In both bending and extension, 
where mlcrostructural damage Induces both damping 
and stiffness loss. The equations of motion are 
formulated using the theory of Internal state 
variables (ISV‘s) to account for damage [10- 
12]. The resulting equations are cast In such a 
way that experimental data can be readily 
Implemented to the model. 

DEVELOPMENT OF THE EQUATIONS OF MOTION 

Consider a free body diagram of a laminated 
composite beam, as shown In Fig. 1. 
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Fig. 1. Free Body Diagram of Lamlnalad Composlle Burn 


The governing equations are 
1) resultant loads [13) 
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(1c) 

2) equilibrium 113], neglecting rotatory 
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3) kinematic [13,15] 


3 U 

c * — 

xx ax 

(3a) 

, 8w o 

e y *~ 

(3b) 

u ■ u 0 + i y z, u « u Q (x), e y * B y (x) 

(3c) 

where equation (3c) is the Euler-Bernoul 11 
assumption. 

4) thermomechanical constitution [10,11] 


e xx * E < e x'x - «1 i T - T r 

(4a) 

q 1 ’ *( k ij * k D1j^ T ,j 

(4b) 


# 1 * “l^xx* *xx* T * *1* 

k 01J ‘ k 01J^°P 
where Is an Internal state variable 

representing matrix cracking In each ply [9]. It 
can be shown that oj Is a kinematic description 

of the locally averaged mlcrocrack geometry 
[10,11]; that Is, 

O, - 4- / V>s (5) 

1 V L S c . C 

where V» Is a local characteristic volume which 
Is large enough to encompass statist Ically 
homogeneous damage. Furthermore, u c Is the 
normal component of crack opening displacement, 
and S c Is the surface area of cracks in V L , as 
shown In Fig. 2. V, may typically be taken as 
one square inch In the x-y plane for laminated 
polymeric composites |3]. For matrix cracking 
the 2 coordinate Is taken to be one ply thickness 
|11]. Equation (4c) Is termed the Internal state 
variable growth law. 


k z 



Fig. 2. Kinematics of Matrix Cracking 


5) thermodynamics | 16,17] 

“l - pc vii- + > r " 0 (6) 

which Is a statement of conservation of energy 
for the material described by equat1ons c (4). 
Although specific forms may be given for uj“ as a 
function of c xx and « 1# this Is covered In a 

later section. Note that although the mechanical 
formulation is one-dimensional in the spatial 
coordinate x, multi-dimensionality must be 
retained in thermodynamic equation (6). 

Now substitute (3c) into (3a) to obtain 


c * — - +2 — ■ 
XX 3 X 3X 
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(7) 




Substituting the above result Into (4a) gives 


>u 0 aa 

a ■ E( + 2 — ^ - a, - ael 

xx v ax ax 1 ' 


( 8 ) 


thickness. This Is generally the case In bending 
and with asyroetrlc layups |15). 

Now substitute (13a) Into (2a) to obtain 


Now substitute (8) Into (la) and (lc) to obtain 

p - £ l* TT * £ i TJ 1 ,/ I 


■ ,h - 


►, * h<*‘ TT> 


(14) 


■ X. 1 E i°n A i ■ 


(9a) 


which Is the governing equation of motion for 
axial displacement. 



, E i°n z i A i ■ 


(9b) 


where “n 1s assumed to be independent of y and 
2 in the 1th ply 111], and 

A* = J f- dA (10a) 

A L 1 

I* s / r z2<1A (1°«) 

yy A L 1 

Furthermore, 

P T s / EoedA (11a) 

A 

H / Eo62dA (lib) 

y A 

The term containing c,. in equation (9b) does not 
go to 2ero even A1 for a symmetric layup 
because is not symmetric (due to the crack 

opening displacements in equation (5)) during 
bending. 


Now, suppose the modulus weighted centroidal 
axis is utilized, such that, 

2*- / S- 2dA - 0 (12) 

A *4 • 

Note that _^ince E Is the initial undamaged 
modulus, 2 represents the Initial modulus 
weighted centroidal axis, which is not damage 
dependent. Therefore, using equations (3b) and 
(12), equations (9a) and (9b) simplify to the 
following for symmetric layups 


£ 1 A 


r au o n 


“u A i - p 


« -li* 

M y * £ i ! yy ^ * ^ E i“ii z i A i * K y 


(134) . 

(13b) 


SlmlUrly. substitute (13b) Into (2c) to 
obtain 


i 2 w 0 n 


r-Z , E i“n' z i A i - H y } ‘ v r (15) 


ix lu i‘yy sx 2 


i-i 


Substituting (IS) Into (2b) thus gives 


aw. 


- p 2 - ir (eA «*> 


(16) 


which is the governing equation of motion for 
transverse displacement. 


Now consider a simple example for the 
constitutive behavior. Suppose that the stress- 
strain relation can be modelled by a Voigt 
element, shown in Fig. 3. 

The governing differential equation for the 
analog is 


xx 


nc 


xx 


Eo8 


(17) 


Suppose we let 
"1 -A 

“11 “ ‘ E i c xxi 

Then equation (17) may be written 

O 1 E(c - Oij • #fi) 
xx v xx 11 ' 


(IS) 


(19) 


which Is equivalent to (44), 4nd (18) replaces 
(4c); that is. 


“11 (t^ * f D l(‘xxi’ ‘xxl* T * °11 )dt * 



•A 

c xxi 


<h> 


( 20 ) 


The term containing in equation (13b) 
represents a shift in thr centroidal axis due to 
matrix cracking, which is nonzero when the matrix 
crack displacements are not symmetric through the 


Substituting (7) and (3b) into (18) gives 


a 11 




5 1 6 X 



0 


stsx 


( 21 ) 
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cr 


FJp. 3. Volpl Model 


T 

£ 


Therefore, substituting (21) Into (14) gives 
s ’ * ™ n 3u 2 ; 3 w 

» I(E i A tt ) + X n Ai A i> irfr ^ ,i Bi 5 i A rr £ 2)J 

W 3t3X 


5U. 


- -£) - P + ±P- 

•t' 3t * r X 3X 


(22) 


C(t) 


E o* £ l 


(t) 


(25) 


where E 0 • E 0 (*j), n l " n l * 

In general, the constitutive »ode1 
represented by equation (4a) Is not In agreement 
with equation (24). However, it Is our Intent to 
show that this assumption Is at least 
approximately correct. *Io do this, first recall 
that equation (24) may be represented by an 
Infinite order differential equation: 


E(t xx - M) * ^ 


A. 


dt 


(26) 


Equation (26) may be written In the same form as 
equation (17), where In this case 




c xx* 


(27) 


Therefore, equations (22) and (23) may still 
be used except that the damping coefficients are 
not constant. For near-constant amplitude load 
histories. It Is now assumed that 


‘xx * *xx(“' C A } 
c xx ‘ ‘xx ( “* e A } 


Furthermore, substituting (21) Into (16) gives 


.2 


* 


.2 n 

0, . t 


. 2 ., 


2^1*yy ?) + 2^ I ) rrr 

3 * 1 3x" 3x" T-l 1 1 A1 


- 7 

+ r, B1 2 . A. — — 2 ) - — ( P A — ) 

e kc X 
-2yT 
-p 

2 ax" 


(23) 


Therefore, if and n fi . are constants In 
time t, the above reduces to linear theory. 

Experimental evidence indicates that 
polymeric composites may not be represented by 
simple mechanical analogs such as the Voigt model 
118]. For small displacements a single integral 
model of the following form has been "shown to be 
accurate [19): 



(24) 


where the kernel function C is a power law in 
time, and is also a weak function of damaae (for 
near-constant amplitude vibrations); that is, 


(26) 


Furthermore, experimental evidence 120) indicates 
that for constant amplitude cyclic loading 
n- 


‘1 * ko A 


(29) 


where k and n^ are material constants* 
For n-»l, the above may be integrated in time to 
give " 

* ^(K, c^) ( 30 ) 

so that equation (27) may be simplified to 

n * n(c^, N) (3^) 

How, consider a more general representation 
of equation (21) given by 


.V 


J ii 




r... 2 

o_ At o 

“ . c te X 


^0 

v Ai 3 x 

^El- ;3w o 

£ 1 1 3t3X 2 


where 

^Ai = r 'Ai (°a , w,N ) * r, Bi * n Bi^ c A ,w,K ^ 
v Ai = u Ai (c A’ u,K )’ v Bi = 


(32) 


(23) 


Substituting (32) into (14) gives 
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aV 


{ X I(1-; A0 )E/^1^(V‘ Itaf) 

' ‘ 7x^Bl he l A ~T * “ax^Bl** „ V 


ax 


3tJX fc 


The last tem In each of equations (34) and 
(37) will produce 2ener-11ke damping due to 
thervonechanical coupling In equation (6). To 
consider these two terms, the fully coupled field 
problem must be solved. That Is, equations (6), 
(34), and (37) comprise a set of three coupled 
equations In three unknowns, u Q , * o , and T. 


7t<’ A TT> - ->x * Sr 


(34) 


where 


P A0 * 


P B1 5 


{.i Va ^ A i 

E / 

n 

I.J w B1 £ 1 A 1 2 1 


hEjA 


and 


n A0 




(35a) 


OSb) 


(36a) 


l . "B1 z 1 A 1 


n Bl 1 


T*1 


hA 


(36b) 


Substituting (32) Into (16) gives 


,' x 2 1(1 * p B2 )E l ! yy * ^2 ^Al h£ i A ’ TT^ 


a 2 - 

> 2 u 0 - 2 

sx 2 ^ r, Al £ l hA 

0\ . C 

sw 

.V 

- fr* 7T> - 



S^W 


siax 


where 


P A1 


M B2 = 


^. 1 w ai e i a i 2 i 

hE^’ 

i.j P B1 £ 1 A 1 2 1 


and 


n Al ; 


E2 1 


W n ^1 ~i A l z 1 
hE^' 

n -2 
Z ul T 'B1 A 1 2 i 


(37) 


(38a) 


(38b) 


(39a) 


(39b) 


yy 


Determination of Material Damping Constants 

The energy dissipation due to damage must be 
at least first order In damage [10); that Is, 

U L " u L (t xx* T * °1 } ‘ u Ll (t xx* T) ‘ °1 (40) 

Therefore, assuming that 2ener damping may 
be neglected, a first order approximation Is 
given by [10,11] 


“L- E ‘xx ; i 

Thus, the total energy dissipation during a 
simple cycle In a specimen of length L is given 
by 


t^+Zr/u ^ 

I J I t 


0 A 


C bi 
XX 1 


dAdxdt 


(42) 


Substituting equation (3) for the case of 
pure bending thus gives 


V 2,/ “ L n 


a 2 v 

/ ; i t^i t 

t, 0 1-1 ’ 1 1 ax 2 


’ll 


dxdt 


(43) 


Thus, substituting (32) into the above 
results in 


/ 

h 


t^+2i/u L 


Li 


. 2 w fl 

£ iVi 


"Bi 2 i _°2 

«k«A 


_B1 - 
£ i 


1 .*2.2 

t U * X 


)dxdt 


(44) 


The above neglects energy losses due to shearing 
[21], which may not be negligible In laminated 
composites [18). 


Now consider the case where the damage Is 
independent of x and no new damage Is introduced 
during the cycle of Interest. It follows that 
for symmetric vibrations, the first term in the 
above equation Is zero, and equation (44) reduces 
to 

n 




-si 


1-1 


, -2 

'BiVi 


(45) 


where 


S = 


t 

f 

t. 


^2i/w 


L 


/ 

0 


s 2 w a 4 -. 


ax 2 at^ax 2 


dxdt 


(46) 


Substituting equation (39b) and rearranging thus 
results In 


‘82 


U 


c 

L 


yy 


s 


(47) 
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Therefore, for * bees of given geometry end 
boundary conditions, equation (47) can be used to 
obtain n B2 * Equation (39b) can then be 

utilized to extract n B1 for certain simplified 

stacking sequences. This Information may then be 
utilized _to predict the damping coefficients 
n B 4 and r> D , as a function of damage In other 

D I D4 

stacking sequences by using equation (39b) for 
other conditions. 

Although It Is In principal possible to 
obtain the’ other coefficients described by 
equations (35), (36), (38) and (39), It Is 

difficult to accurately measure u£ during 

dynamic axial tests, so that and ere 

not easily evaluated by the method described 
here. However, the_ stlf fness loss coefficients 
w A0 , Z B y v A1 . *nd Z B2 may be obtained from 

direct static testing |11,12). 


Experimental Determination of Damping 


controlled laboratory experiments, the relative 
change In damping due to a progression of 
mlcrostructural damage *111 be relatively 
Independent of the system damping which should 
remain constant. 



Because of its relative simplicity, a 
cantilevered bear, has been selected as the 
current test configuration for studying material 
damping. (Jt should be noted that the 
theoretical development described In the previous 
section Is valid for all beam boundary 

conditions. We are also Investigating other bear, 
boundary conditions such as the free-free beam 
but will confine the present discussion to the 
behavior of a cantilever beam). The laboratory 

setup is shown schematically in Flo. 4. A coupon 
1.0 in. (25.4 rm) wide and 12.0 in. (304. B mm) 
long Is clamped In a beam bending fixture to 
achieve a cantilever of the desired length. The 
beam, is set into vibratory motion by a motorized 
wneel which rotates until a post engages the free 
end of the beam,, achieves a preselected 
oeflection and then releases the free end. 
Therefore, the frequency and amplitude of the 
beam are controlled by varying the length of the 
beam, and the initial deflection. The free decay* 

response Is recorded by a strain gage mounted on 

the beam near the clamp. The bend fixture is 
housed in a vacuum chamber and the test is 
Initiated externally by a computer which also 
acquires the strain oace data in digital form. 
The output of the strain cage is recorded at 1 50*3 
data points per second. The amplified signal 
from the ' 1000 ohm strain gage is filtered to 
Improve the data reduction accuracy. The damping 
of the beam Is calculated by the logarithmic 
decrement method. 

In the present study the outer fiber strain 
amplitude was selected so that the logarithmic 
decement response of the beam Is linear. The 
system, (joint) damping has been estimated from 
tests on several well known materials and has 
been subtracted from the experimental damping 
values measured from the composite. It is 
recognized by the writers that the cantilever 
beam may not be the most desirable test 
configuration for a study of material damping. 
However, our primary objective is to study 
camping as a function of mlcrostructural 

damage. Therefore, by constructing carefully 


Ftp, 4. Schematic Drawing el Tailing Apparatus 


Analysis of t Cantilevered Beam 

This section provides the specialize! ion of 
the damage-dependent beam theory for the 
cantilever beam.. The analytical procedure for 
ostermininp the materiel damping constants from, 
the experimental data is also described. 

For an undamped prismatic cantilever beam, 
undercoinc free vibrations, the solution is of 
the form. 1 22] 


v 

o 


(t.X) 


y 

r.=l 




where 

f r (t) = C, sin ♦ CnCOS * f 


and 


(48) 


(49) 


cos a 1 + cosh a l 

;; i E 

m l v sin £ L + s inn a 
m m 

(sinh a x - sin a x) 
m rr. ' 

+ cosh ax - cos a x] 
jr m 

AT so 

a L - 1.875 rr— 1 

(m-l/2)r m>i 


%.(*) 


(50) 


(SI) 


259 


ORIGINAL PAGE IS 
OF POOR QUALITY 







and 

- (0.597,/L) 2 |(]-p B2 )E 1 I* y /pA) 1/2 m-1 

- I (m-l s ),/L) 2 [(l-w B2 )E 1 I^ y /pAl 1/2 m>l (52) 


deflection an/tj and C? depend on the Initial 
conditions. Equation (52) may be used to 
determine _ the laminate bending 

coefficient Z^ when w n can ^e obtained from an 

experimental test. 

In our Initial phase of research, the 
damage-dependent constants are more easily 

determined If we consider only the first mode of 
vibration and confine our experiments to the 
first mode. (Later research will Include other 
modes of vibration.) Considering only the first 
mode and a beam with the following Initial 
conditions: 


- / 
0 


2, - /u l f L 2 ,1^875/ 


; - 


5.2734 


2 

2 cos wjt 


(0.7341 (slnh 1.875 f + sin 1.875 £ ) 


+ cosh 1.875 i «■ cos 1.875 *-] 2 dxdt (58) 

Integrating the above equation and 
substituting Into equation (47) thus gives 


n B2 ’ 


“1* 


1 


(59) 


20 ' 121 yy L a 

The energy dissipated per cycle, u£, of a 

vibrating beam Is related to the experimental ly 
determined logarithmic decrement, a, by the 
following relationship 


w e (0.L) - t 

(534) 

l « 

1 

2 U $ 


(60) 

sw 

jt (°* x > “ 0 

(53b) 

where 

beam. 

u s 

The 

Is the maximum 
strain energy of 

strain energy of the 
the beam Is* given by 


Then « 0 and C 2 * 1- Furthermore, 

♦ : (x) • A : l (0.7341) (sinh 1.875 £ - sin 1.875 *) 

+ cosh 1.S75 i - cos 1.S75 i| (54) 

Substituting the above into (48) gives 
* 0 (t,x) - AjCos I (0.7341) (slnh 1.875 £ 

- sin 1.875 ^ + cosh 1.875 £ 


U « i J L J E (-L^-) 2 2 2 dAdx (61) 

A SX 

and for a cantilever vibrating in the first fr -vJe 
with an initial tip deflection, 4 , 

U_ - 2.2016 -4- l E^.T 2 (62) ' 

5 l i=l 1 1 

Substituting equations (108), (50), and (62) Into 
(59) gives 

"62 " °- 31S27 (—“J 


- cos 1.875 fl (55) 

Satisfying Initial condition (£3) gives 


'1 5.2734 


(56) 


Therefore, 

w 0 (t.x) - 5 ~ 273 ^ COS [0.7341 (sinh 1.875 f 
- sin 1.875 i) + cos 1.875 j- 
- cos 1.875 j|l (57) 

Substituting (57) into (46) thus gives 


The procedure used herein to determine the 
damping constants, t^, is to experimentally 

measure the logarithmic decrement for a 

[0/90/0] laminate and calculate n B1 from 

equations 5 (63) and (39B). Since we ire only 

considering matrix crack damage at present, we 
will assume that all the damping may be 

attributed to the 90* plies of the [ 0/90/0] s 

laminate. This will provide an experimental 
measure of the damping in a single 90* ply with 
matrix crack damage. Substituting equation (39B) 
Into equation (63) and rearranging gives the 
following expression for the logarithmic 

decrement 


3.142 u 


n 


T 2 


l. r, Bi Vi 


1 


y f.a>2 


- 2 


(64) 


i = l 


i°ri 


Once r IQ< j is determined as a function of damage, 
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eouatlon (64) can be used to predict the damping 
of a cantilever beam of any cross-ply laminate 
stacking sequence. While equation (64) appears 
linear In damping this Is not the case because 
the damping constants, are also function of 

frequency. (See equation 33). 


Comparison of Experimental and Analytical Results 

Experimental results have been obtained from 
cross-ply laminates with various combinations of 
0* and 90* plies of A5-4/3502 graphite/epoxy. 
The lamlnetes were fabricated using a standard 
prepreg tape layup and hot press curing process 
following the curing procedure recorrmended by the 
prepreg tape vendor. The approximate per ply 
thickness of the cured laminates was 0.005*' 
(0.127 mm) and the- fiber volume fraction was 
approximately 62Z. The measured lamina 

properties are given In Table 1. These 
properties are typical of this material system. 

The theory presented herein was developed 
only for material damping. Therefore, it was 
essential to separate the system damping from the 
measured values of logarithmic decrement to 
determine material damping constants and to 
compare theoretical predictions to experimental 
values. Using test results for a variety of 

materials, the system damping for the current 
configuration of a cantilevered beam vibrating 
predominantly In the first rede has been 

determined to be approximately 3.52x10"’. This 


Table 1. 

Material Properties for Hercules 
A54/2502 

Lamina Prose 

rtles 

r 

-1 

21.0 x 1C 6 i 2. OS psl (144.6 GPi) 

r 

w 22 

1.2? x 10 6 i 2. IS psl (?.5S EPi) 

1 

b. - * 
iC 

0.694 x 10 6 psi (<.7? E?4) 

v :2 

0.210 ± 2.7S 

F tu 

326000 i 3.5S psl (2.24S KPj) 

r tu 

11025 ± S.ES psl (76.4 HPi) 

c tu 

0.0144 i 4.6S in/in 

c tu 

0.00772 * 6.7S in/in 


* Theoretically calculated from test data for a 
[±<5)2 S laminate. 


value of system damping was considered to be a 
constant for all laminates in this study. 

Considerable attention was given to the 
experimental test parameters that yielded the 
most repeatable results and held constant those 
parameters which influenced damping such as 
frequency and stress -amplitude. This procedcure 
allows for the most direct comparison between 
theory and experiment for the laminates discussed 


herein. All test results reported herein were 
obtained »t the freouency of the first »ode of 
vlbntlon. The Initial tip deflection we s 

selected so that the outer fiber stress near the 
root was within the stress amplitude range where 
denying was constant and the free decay was 
linear when plotted on a log sc&le, as shown In 
Fig. 5. 




8 t s 


- [o/*D/§ Cw— »*«Ax"Tr 


ac *c tc ioc 


Vm ttocfj 


Since the contribution to material damping 
of the matrix material Is considered to be 

substantially greater than that of the graphite 
fibers, as a first approximation, the assumption 
Is made that all the material damping of the 
cross-ply laminates may be attributed to the 9^ 
layers. Furthermore, we are confining ^the 
present studv to the single damage node of matrix 
cracks in the 90* plies. Tnerefore, it is 

necessary to specify only one material damping 
constant, r>^, in eouation 39n, to determine the 
damping of D any cross-ply laminate. Using the 
experimental results for the (0/90/0) c lamina. e, 
r, B1 for a single 90* ply was calculated from 

equations (63) and (39b) to be -.86 Ibs- 

sec/1n^ for an undamaged laminate. This value 

of r^. and equation (64) was used to predict the 
undam.coed material damping of the 1 0/93 J • 

I0/9D-L and l0/90 2 ] 3 ^ laminates. The comparison 
of the predicted^ values of damping to the 
experimentally measured values for the undamaged 
laminates is given in Table 2. The predicted 
values are within ±SS of the experimental 
values. This is consistent with the range o« 
repeatability of the experimental data as 
measured by the standard deviation of numerous 
replicate tests. 

A complete characterization of the damage 
dependent material damping of oraphi te/epoxy 
laminates is in progress. The experimental 
results confirm a strong relationship between 
material dampina and microstructural carnage. 
Figure 6 shows the relationship between ^material 
damping and microstructural aamage- . r ' cure 6 
shows the relationship between material damping 
and matrix crack carnage in the 90* layer of a 
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Table 2. 

Predicted and Experimental Values of 
Oamping Determined by the 
LogarltNnlc Decrement Method 

LAMINATE 

FIRST 

EXPERIMENTAL 

PREDICTED 

STACKING 

SEQUENCE 

MODE 

FREQUENCY 

DAMPING 

DAMPING 

|0/901 2j 

18.5 

4.52xl0* 3 

4.20xl0‘ 3 

|0/90 3 ] s 

IB. 7 

5.1U10* 3 

5.56xl0* 3 

I0/90 2 ] 3s 

25.3 

10.64xl0* 3 

10.«xl0* 3 

10/90 3 ] s 

laminate. 

(The damage 

states were 


developed In uniaxial tension-tension fatigue.) 
As can be seen In Fig. 6, damping, normalized by 
the Initial value of damping for the undamaged 
laminate. Increases by about a factor of 2 as the 
matrix crack damage state progresses to the crack 
saturation state. 


Summary and Conclusions 

Equations of motion for a vibrating beam 
have been developed for composite materials which 
account for the Influence of microstructural 
damage on beam stiffness and material damping. 
This new theoretical approach introduces damage- 
dependence into the constitutive relationships 
via Internal state variables which characterize 
the current microstructural damage state. Since 
the constitutive relationship Is time dependent, 
the equations of motion contain explicit terms 
representing the damage-dependent material 
damping. The equations of motion are developed 
for a laminated composite beam by employing 
standard lamination theory and assuming that 
Individual ply contributions to material damping 
can be represented by a rule-of-mixtures 
formulation. The equations of motion may be 

solved to determine the damage-dependent natural 
frequencies and mode shapes of the beam. The 
beam theory is also used to construct a 
mathematical model of damage-dependent material 
damping at the laminate level which Is the 
primary subject of this paper. 


damping can be attributed to the 90* layers. The 





undamaged damping constants for a 90* ply were 
determined from the experimental data of the 
10/90/0] s laminate. The undamaged damping of the 
I0/90U* |0/90 3 ] s and (0/90 2 ] 3s laminates was 
predicted by the theoretical model of damping. 
The predicted values for these laminates were 
within ±8* of the measured values. While damage- 
dependent data are not fully developed at this 
time, the excellent agreement between the theory 
and experiment for the undamaged damping provides 
verification of the theoretical formulation. 
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8.3 An Experimental Investigation of Damage -Dependent 
Material Damping in Laminated Composites 
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DAMAGE-DEPENDENT MATERIAL DAMPLNG LN 
LAMINATED COMPOSITES 


S. A. Smith • 
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C. E Harris 
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Abstract 

An experimental program is reported in which the 
material damping of laminated graphite /epoxy specimens 
was measured as a function of matrix cracking. The damping 
test apparatus is designed to measure the first mode free 
vibration response of a cantilevered beam. Damping is 
calculated by the logarithmic decrement method. Matrix 
cracks are introduced into the material, and damping is 
measured at several levels of damage. The results from six 
dif/ereni laminate stacking sequences are recorded and 
compared to each other. Damping is found to be significantly 
*Uecicd by the laminate ply sucking sequence, and is also 
seen to increase dramatically with the presence of matrix 
cracks. Finally, mathematical model predictions are 
compared to the experimental results. 

Nomenclature 


A 

• beam cross -sectional area 


- cross-sectional area of the ith lamina 

£ 

- axial modulus of elasticity 


- reference axial modulus of elasticity 


- axial modulus of elasticity of ith lamina 

3 yy 

- modulus weighted moment of inertia 
about y axis 

My 

- beam moment resultant about y axis due to 
temperature change 

N 

- number of cycles of applied load 

r 

- number of plies in laminate 

x>z 

- beam transverse load per unit length 

S 

- term used to obtain laminate damping 
coefficient 


- temperature 

t 

- time 

Ui 

- total interna] energy’ dissipation due to 
matrix cracking 

Us 

* maximum strain energy of a deflected 
beam 

w O 

- transverse displacement of neutral surface 

X 

- beam axial coordinate direction 

= 

- beam transverse coordinate direction 

z i 

- o component of distance from centroid of 
ith ply to modulus weighted centroid 

C 

- coefficient of thermal expansion 

C} 

- axial component of matrix cracking 
damage internal st 2 te variable (ISV) 

c li 

- axial component of matrix cracking 
damage ISY in ith lamina 

£ 

- beam tin defection 
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b 

- logarithmic decrement 

‘xx 

- axial strain component 


- damping coefficient 

T 1B» 

• bending damping coefficient in ith lamina 

MBi 

- bending stiffness loss coefficient in ilh 
lamina 

*3 

- modal frequency of the 1st resonant mode 

P 

• mass density 

°A 

- stress amplitude of cyclic loading 


Introduction 

Microstru crura) damage in laminated fibrous composite 
materials has been found to alter several material 
properties^. * Extensive modeling and testing has been 
performed to determine stiffness loss in laminated composites 
as a result of matrix cracking and interply celamin 2 t>ons. In 
addition, several investigators have noted and measured the 
changes in camping caused by carnage 4 "^. However, no 
attempt has been made to quantify or mode) the effect of 
carnage on material damping. Tne objective of this study was 
to define a direct relationship between microstructural 
damage in a laminate and the resulting change in the damping 
ca parin* of that laminate. 

Tnis paper presents and discusses the results of art 
experimental program which has investigated the effects of 
matrix cracking on material damping of laminated 
graohite/epoxy specimens. Tne effect of ply stacking sequence 
on material damping is also discussed. A mathematical model 
developed by Alien, Karris, and Kighsmith' has been 
applied, and the resulting predictions of damage-dependent 
damping are compared to the experimental results. 

Mathematical Model 

Many mathematical models have been developed which 
attempt to predict the materia) damping of composite 
materials 1 However, none of these models are able to 

predict the material damping of composites as a function of 
mi crostru crura) damage. Tne model discussed in this section 
was developed by Allen, Harris, and High smith' at the 
Mechanics and Materials Center at Texas A£cM University. 
Tnis section provides only a brief overview of this model, as 
detailed derivations and explanations appear in Reference 9. 

Tne model uses the concept of continuum carnage 
mechanics, using internal state variables to account for damage 
in the constitutive equations. Tne internal state variables 
may be defined for each particular form of damage in the 
composite. Tne constimtive equation for the material system 
used in this stutiv is based on a nonlinear Voir, model 
representation, shown in Figure 1. Tne general consecutive 
equation for this model is: 




--CAT) 
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where ti is a conaunt material property, it the coefficient of 
thermal expansion, and AT is the change in temperature. 
Experimental evidence has shown that polymeric composites 
are not accurately represented by this simple model. In the 
model development, it is shown that a Voigt model in which tj 
is not assumed to be constant for all loadings, but is instead a 
function of initial stress, c A , frequency, «, and the number of 
cycles, N, will closely approximate Ihe actual material 
constitution of polymeric composites for small displacements. 



Figure 1. Voight Model 

The general thermomechanical constitutive relationship 
for a material, including the internal state variable C], is : 


Cjq. * HLjq. - a] - c AT) . (2) 

Equating (1) and C2), the internal state variable 

a : - Jt* (3) 


for a single ply within the laminate is determined. 

The equation of motion for transverse displacement 
.(bending) is developed using the constitutive equations, 
equilibrium, and strain displacement relationships for an 
Euler-Bemoulli beam. 


J’vy -j • Ej a li H ) 1 

i=l 


o 

c-w fi 


n 

V 1 


3, * **'<>* 


«> 


where 


and 


l 'yy“ J^ dA 

A 





a (AT) z dA 


(4a) 


(4b) 


Replacing t xx in (3) with displacement variables, and 
introducing > coefficient Pbi ioT > ,h PV- which de * cribri 
the loss in (tiff ness caused by matrix cracking. Ojj for a pure 
bending case b ecomes: 

ch~ D Tibi 

•u-^fcT' T~ 



The model development also accounts ter axial displacement, 
but will not be considered in this study since only the bending 
cue has been investigated experimentally. 

Substituting (5) into (4) yields: 


£ 
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where 


and 



]^UBi Ej V 


MB2 


i-1 


El i 


yy 




tib2 


i*l 


lvv 


( 6 ) 


(6a) 


(6b) 


The determination of the laminate damping coefficients 
is accomplished by calculating the energy dissipated per cyde 

of vibration, U£. Assuming that Zener and shear camping 

effects can be neglected, that the damage distribution is 
uniform over the length of the specimen, and that additional 
damage is not introduced during the cyde of interest, the 

bending damping constant, fj£ 2 , is approximated by 


where 



n 

2/lbi n 

i*l 


(7) 


(7a) 


and 


f r g»o 

J J ers?- di~c 


^rdxdt 


(7b) 


For a cantilevered beam vibrating in the first mode, w 0 is 
derived as 


Wot,*) = M 0 - 7341(sir ' h U75 f - *** 1 - S 75 r ) + 

cosh 1X75^- cos 3£73^ l6) 

where the constants were determined by the boundary’ 
conditions. Substituting (S) into (/d) and integrating. 
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Tabic 1. Lamina Properties for AS4/3502 Graphite/Epoxy 


S =20.12-7^- . <*> 

where & is the initial beam tip deflection and L is the length 
of the beam. Substituting (9) into (7) yields the following 
expression: 

6>i 

l£- 20.12 T lB2lyy"]T" - (10) 

The logarithmic decrement. « defined to be one-half of the 
energy dissipated per cycle, , divided by the maximum 
strain energy of the deflected beam, Uj. For the given 
configuration, 

n 

U, - 3.2016 ^Ej A; . (11) 

L »«1 

Therefore, 

n 

5k3.K2(di ————— , ( 12) 

n 

X £ i A i v 

i«l 

• which is a first order approximation of the logarithmic 
decrement for a symmetric, cantilevered laminate, vibrating in 
the first mode. Additionally, this form assumes that the 
displacements are small, that the beam is not shear 
deformable, and that lamina property' 7^5, is a function of c^, 
to, N, and ply orientations. 

Tne values of tjbj for a pven damage state, initial stress 
level, and frequency are determined by measuring the 
logarithmic decrement and extracting 7^ for each ply. It is 
also assumed in this model that, for cross-ply laminates, 71 gj 
for 0 degree plies does not change when matrix cracks are 
introduced in the 9D degree plies. Therefore, for this study, 
TjBiCO) * De assumed constant for all damage levels. 
However, Tigj^Q) be defined as a function of matrix 
cracking, and will be determined experimentally. 

Experimental Program 

Specimen Preparation 

In order to meet the objectives of this study-to isolate 
matrix cracking as the only damage mode-cxoss-ply laminates 
were chosen to be tested. 12 inch square plates were prepared 
from Hercules AS4/3502 Graphite/Epoxy' preimpregnated 
tape, and 1 inch wide specimens were subsequently cut from 
these plates using a high speed diamond -edge cutting wheel. 
Table 1 lists the ply properties of the material which have 
been obtained from appropriate testing methods. Table 2 lists 
the laminates selected for testing in this study. 

To minimize undesirable damage caused by the grips of 
the loading machine, a thin coating of epoxy was spread over 
the ends of the specimens . Tne specimen edges were then 
polished to insure high quality edge replicas. lOODD single 
ejement strain gages, aligned with the beam axis, were bonded 
to the specimens to measure the response signal during the 
camping measurements. Tne final specimen configuration is 
seen in Figure 2. 


Prppmy 

y*3vr 

til 

2 U»* lO^psi 

*22 

139 x lO^psi 


0.694 x 10 6 p® 


0310 

fiber Volume Fraction* 

62% 

Calculated from test data for a I ±4 5} 2 * 

laminate 

Manufacturers Pre-Cure Value 


Table 2. AS4/3502 Laminate Types Tested 

Laminate 

Stacking Sequence 

A 

10/90/0), 

B 

190/0/90), 

C 

IO/ 9 O 3 ), 

D 

I 9 O 3 /O), 

E 

10/90)2, 

F 

150/0/502), 



' Figure 2. Specimen Configuration 

Pamrring- Measurements 

Because of its simplicity, a cantilevered beam 
. configuration was chosen for the camping apparatus. A 
• support was designed which minimized specimen slippage at 
: the boundary, and avoided inducing damage in the specimen. 
Figure 3 is a scale drawing of the fixture. To conduct damping 
' tests, the specimen was placed in the slot with a 
predetermined length extending from the support. This length 
was chosen to produce a desired beam frequency. 

Many experimental techniques have been employed in 
the past to study the material damping of composites. One 
concern in these studies has been the effect on damping caused 
by external factors 5 ' 6 ' 16 "^. Some of the undesirable damping 
contributions encountered in past investigations are 
atmospheric pressure, friction caused by contacting instruments 
- and/or boundaries, and additional masses. Regardless of the 
method used for damping measurement, considerable difficulty' 
has been encountered by almost all experimentalists in 
obtaining repeatable and accurate results. Since this study 
focused on the effect of damage on material damping, all 
other damping contributions were minimized and held constant 
to establish an exact correlation between damping and 
damage. Atmospheric pressure has been found to influence 
camping. To eliminate this effect from measured damping, all 
tests were conducted inside a vacuum chamber. To study the 
effect of the epoxy end tabs on material camping, tests were 
conducted firs: without the end tabs, and then with the end 
tabs on the specimen. Tne measured damping for a given 
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configuration ^ecraucd with the addition of the epoiry md 
UD>. T>us WAS. Attributed to the resulting decrease w frequency 
roused b ) 1 the smdl mass At the tree end. And did not Alien the 
Actual material damping capacity of the speorrenv. 




Tipu X loundcy Swppon 


In order to exdte the beam at a prescribed amplitude, a 
mechanical system was developed to “pluck" the free end of 
the beam to a given deflection. An electric motor is attached 
to an aluminum mount which slides in a groove cut into the 
baseplate of the specimen support. Also attached to the mount 
is a routing wheel which is powered by the motor. An arm 
extends from the wheel and is used to contact the free end of 
the specimen- The excitation source and support are shown in 
Figure 4. Upon initiation of a damping experiment, the wheel 
is routed until the arm contacts the specimen. The arm then 
continues to route, causing the beam to deflect, and eventually 
releases the specimen into a free vibrating motion. Tne tip 
deflection of the beam is adjusted by moving the contact point 
of the arm to the specimen. 



Figure 4.. Damping Excitation Apparatus 

Tne strain response signal was transmitted to a half 
bridge contained in a Micro Measurements EAM-1, which was 
also used to indicate static strain values. In order to obtain an 
initial strain level that produces a signal of sufficient 


Amplitude lor the damping measurements, the arm-to- 
spe o men contact point is moved by adjusting the position of 
the wheel mount on the baseplate. Figure 5 is a schematic of 
the test svstem. The entire process of specimen cx citation and 
data collection is controlled by a personal computer. An analog 
and digital interface board was installed into the computer. 
The control of the motor and the collection of the dau is 
provided by software. The daU sampling rate was adjusted for 
each lest to collect exactly 50 dau points per cycle of beam 
vibration. 


|w» » /— * 



To accurately measure the damping and frequency of the 
specimen, a sufficient number of cycles had to be captured. 
Since dau was collected at 50 points per cycle, capturing 9000 
dau points results in recording 180 cycles. This number of 
cycles was adequate to describe the decay of the specimen 
; response in vacuum. With the number of data points and the 
sampling rate prescribed, the computer was prepared to start 
the damping experiment. An initial zero load strain reading 
was made which was subtracted from the output signal. After 
this value was recorded, a 5 volt analog signal, generated by 
the D/A interface board, was supplied to an external circuit 
board. This circuit then directed a 12 volt power supply to 
transmit a signal to the motor, causing the mechanical arm to 
rotate. At this point, the computer began to read the dau 
signal from the stationery beam. As the rotating arm 
contacted the specimen and deflected downward, the resulting 
tensile strain was read by the computer as a positive voltage. 
When the arm released the specimen free end and the beam 
deflected upward, the strain approached the initial unloaded 
value. Once the beam crossed this initial deflection plane, the 
strain became compressive, and was read as a negative value. 
As soon as the computer recognized the negative value, it 
ceased to supply a voluge to the external circuit, causing the 
wheel to stop routing. Then the computer collected dau at 
the prescribed sampling rate. By collecting data only after 
the motor had been turned off, the electrical signal and 
vibration caused by the motor, would not influence the data. 

One measure of the logarithmic decrement, 5, is given by 

5 ..In (ill) , (13) 

x i 

.Where xj is the peak strain amplitude of the cycle. For 
most materials, however, the peak decay of consecutive cycles 
is extremely small, and difficult to resolve from experimental 
data. Another approach is to measure the logarithmic 
decrement between a finite number of cycles. Tne logarithmic 
decrement becomes: 

. 04) 

However, if the N“ n dau point is not a valid representation 
of the response at that particular cycle, error in c can be very' 
significant. If the decaying response is not greatiy influenced 
by the presence of shear deformation, friction, or atmospheric 
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pressure Us discussed earlier), the free decay will be 
exponential, and independent of strain levels. By determining 
these effects to be negligible, the logarithmic decrement was 
measured by fitting a straight line through the natural log of 
the peak values. This was done for the positive (tension) and 
negative (compression) peak values. Figure 6 is the positive 
peak value decay for a 10/90/ 0) s specimen plotted on a semi- 
log scale. The degree of "linearity", and the samll amount of 
experimental data scatter of the decay can easily be seen in 
this gTaph. 



Cydc Number 

Fifure i. Fiu V«3t*« Cwciyior Cr/lp 


Damage Inducement . 

The objectives in the damage inducement phase of this 
study were to introduce matrix cracking in the specimens, and 
to quantify the resulting damage. A 20 kip MTS £10 testing 
machine was used to load the specimens. All of the specimens 
were loaded at a rate of 20 pounds per second as recommended 
by ASTM Standard 3035-76. A uniaxial extensometer which 
was attached to the specimen with Tesaiix double-stick tape 
was used to measure strain in the specimens. 

Edge replication was used to document matrix cracking. 
Since matrix cracks generally extend across the width of the 
specimen, as indicated by x-ray radiographs, visible cracky 
can be detected on both sides of the specimen. Edge replicas 
indicate the location of cracks with respect to the thickness 
. dimension, information which is not easy to obtain using 
radiography. However, edge replicas do not reveal matrix 
cracks in the 0 degree plies (axial splits), or areas of 
aelamination. Previous studies of ASi/2502 cross-ply 
laminates show a distinct damage progre s sion for monotonic 
loadings, in which matrix cracks form only in the 90 degree 
plies at relatively low loads. As the load is increased, the 
number of matrix cracks in the 90 degree plies grows 
substantially, until a matrix crack saturation level is reached. 
Higher load levels will not produce additional cracks in the 
90 degree plies. However, axial splits and interlaminar 
delaminations may form, depending upon the cross-ply 
stacking sequence. Since this study was focused on matrix 
cracking, all tests were conducted below the matrix crack 
saturation level, determined for each laminate, minimizing 
the presence of higher damage modes. 

Several measurements of the stiffness of the undamaged 
specimen were made by loading the specimen to the 500 pound 
level, and recording the resulting stress and strain. 
Subsequently, the load was slowly increased, and edge 
replicas were taken at several intervals. When matrix cracks 
first became detectable, the load was reduced and the stiffness 


was measured again. Because matrix cracks are not always 
spaced uniformly along the length of the specimen, and since 
spacing may be different in non-consecutivc 90 degree plies, the 
number of cracks per inch in each ply was averaged over a 2 
inch gage length. The page section was the two inches closest 
to the root of the specimen when in the clamped damping 
apparatus. The ply dimensions varied slightly from specimen 
to specimen. Therefore, damage was quantified by crack 
surface area per inch of length. A numerical value was 
determined by multiplying the cracks per inch of a ply by the 
cross-sectional area (thickness x width) of the ply containing 
the crack. The surface area is approximately twice this value 
since each crack has two faces. A diagram of a matrix crack in 
a single ply appears in Figure 7, which indicates the 
dimensions of crack surface area. 
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Ey?ERJMEKTAl RESULTS 

The first step in conducting the damping and damage 
experiments was to measure camping for 10jg and [90jg 
laminates, which would be used to develop the undamaged 
damping constants for 0 and 90 degree plies. Using the average 
damping and frequency values of £ tests, as well as the 
specimen dimensions and properties, the model constants for 0 
and 90 degree plies were calculated. Tne constants are in units 
of lb-sec/in 2 , and are: 

RSi(O) * 164*625 
*l5i(90) * 1*6**4S . 

To illustrate the difference in the damping of these laminates, 
the peak value decays are plotted together in Figure £. 



Cycle Nunber 
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These values w err used it the model to predict the 
undamaged damping values tor each of the specimens tested in 
this studv. Figure V compares the measured and predicted 
damping values for each of the specimens. The measured 
values represent the average of 5 damping measurements taken 
for each specimen. The next step in the experimental 
procedure was to measure the logarithmic decrement for a 
|0/90/0] s specimen as a function of damage in order to develop 

a relationship between T}Bi(90) the matrix crack surface 
area, S. 



ip.oiw a : jc r c oc ds cc ec r. r 

Umiu po/o/k; io/sq; 
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Laminate A. fQ/90/Qlc 

This particular stacking sequence was chosen to develop 
the damage-dependent damping relationship for matrix 
cracking in a single 90 degree ply because the 90 degree ply was 
not adjacent to any other 90 degree plies, thereby restricting 
the si zt of the matrix cracks to one ply thickness. Once a 
relationship was obtained for the damping of the specimen, 
the mathematical model parameter, Tig: was calculated as a 
function of matrix cracking , and subsequently used to predict 
damping for other AS4/3502 Graphite/npoxy cross-ply 
laminates at various levels of matrix crack damage. Figure 10 
is a plot of the measured logarithmic decr em ent vs, crack 
surface area for this laminate. The value Tfg^Q) as a 
function of crack surface area was calculated from the model at 
each damage state for this specimen, and is seen to follow a 
curve of the form y=Ax°. A power law curve fit was applied., 
and is shown in Figure 11. This approximation is valid for 
crack surface area values greater than 0.005 ir.^, which 
corresponds to 1 crack per 12 inches of specimen length. 
Damage below this level was assumed to indicate an 
*■ undamaged" specimen. Therefore,, for the case in which 
matrix cracks could not be detected, was assumed to be 

the value extracted from the damping test of the [90]g 
laminate. This approximation was then used to predict 
damping of all of the cross-ply specimens at any given matrix 
crack damage state, including the undamaged state. 

The reduction in stiffness was extremely small for this 
laminate. Because the data scatter was great in the stiffness 
measurements, and the reduction was minimal, no quantitative 
relationship between axial stiffness and matrix cracking was 
established. Figure 12 is a plot of normalized stiffness as a 
function of matrix crack surface area for this I0/90/0] s 
specimen. Matrix cracking in a laminate is seen to cause only a 


slight decrease the axial modulus (less than 1% for this 
laminate at saturation), whereas the damping increase is very 
wpnificanl (3751 for the Mmt level of damage). Therefore, 
damping cin be considered u» be i superior method of non- 
oestructively evaluating the matrix cracking in a laminate. 
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Figure 11. NormAliuKf AuaJ Stiiinmi v%. 7oa0 Mavis Crack Arm* 

io/»/o: t 


Lemingtg Li£QZCZMs 

The next laminate tested was also a six ply laminate, 
with a single 0 degree ply between twe 90 degree plies in each 
half-plane. Tne undamaged bending stiffness was predicted to 
be significantly less than that of the l0/90/0) s laminate, 
which resulted in a smaller fundamental frequency for the 
same length beam. Therefore, the length of the [90/0/90J s 
beams was shortened to increase the frequency, in an attempt 
to match it with that of the l0/90/0] 5 laminate. However, the 
specimen had to be long enough to minimize the presence of 
shear deformation^. Tne frequency for the undamaged 
!0/90/0] s specimen was 20.7 Hz, and the frequency of the 
]90/0/90] s beam was KJ Hz. _ 

The mathematical formulation for the logarithmic 
decrement in 04) includes a term which defines the damping 
to be a linear function of frequency, for a constant value of Tig;. 
Goser inspection of the model development reveals, however, 
that Tjgj is not constant with frequency, and in fact, the 
logarithmic decrement is predicted to follow a Zener like 
curve with respect to frequency. Tne damping test apparatus 
does not provide for testing a single specimen at a wide range 
of frequencies. Since the differences in frequency for all of the 
tests are not dramatic (less than 8 Hz difference), the 
differences in damping for the various frequencies are 
approximated by the linear relationship of the model. 
Therefore, for this study, rjgj was assumed to be constant with 
respect to frequency. If specimens at much higher frequencies, 
i.e., greater than 50 Hz had been measured, this assumption 
would not necessarily be valid. 

Tne measured and predicted damping values for the 
undamaged I90/0/90] s laminate can be seen in Figure 9. Matrix 
crack initiation was found to occur at approximately 1200 lbs. 
for this stacking sequence. In this laminate the initial damage 
suite was one in which matrix cracks were abundant in the 
outer 90 degree plies, and minimal at the mid-piane at a load 
of 1405 lbs. At a slightly higher load of 1496 lbs., the cracks in 
the inner 90 degree plies increased by a fa cor of 12, while the 
cracks in the outer plies increased by less than 2 factor of 3. At 
the saturation level, all of the 90 degree plies had 


approximately the same crack spacing (about 30 cracks per 
inch). 

The mathematical model predicts that the damping will 
be affected more by matrix cracks far away from the mid- 
plane, than matrix cracks which are near the mid-plane. 
Therefore, it would not be an accurate depiction of damage- 
dependent damping to relate the damping strictly to the total 
crack surface area in a Laminate. The location of the cracks 
with respect to the mid-plane must be taken into account when 
reporting the data. The reported crack surface area was 

therefore weighted by a non-dimensional factor ( Ij/tP, 

where £j is the distance from the centroid of the crack to the 
mid-plane, and t is the laminate thickness. The equation for 
this weighted crack surface area, S*, is 

N 

07) 

where q is the number of cracks per inch of length in the i 1 * 1 
ply, and tj is the thickness of the i lh ply, and w is the width 
of the specimen. The factor of 2 accounts for the number of 
crack faces per matrix crack. It is stressed that this value was 
not used in the model to make damping predictions, but is 
simply used to identify the amount and location of matrix 
cracks in a laminate. Figure 13 indicates the relationship of 
measured and predicted damping to S* for this specimen. 
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■' Laminate C. fO/9CV& L 

This laminate contained 6 consecutive 90 degree plies 
sandwiched by single 0 degree plies on each outer surrace. Tne 
damage progressed steadily from crack initiation at 32 ksi, to 
saturation at 62 ksi. Tne damping values increased 
significantly for this specimen, and the measured and 
predicted logarithmic decrement is plotted against 5* ir. 
Figure K. Although at first glance there appears to be a great 
deal of data scatter in the camping measurements taken at 
the higher damage levels, further inspection of the scaling 
reveals that the greatest amount of variation in camping 
measurements at a single damage state is less than 4 percent. 
Tne model predicts the increase in logarithmic decrement 
accurately for this laminate at all of the damage levels. 
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This laminate scries had the smallest bending stiffness of 
all of the laminates considered. The natural frequency for the 
undamaged specimens was 32.7 Rc, which was also lower than 
that of all of the other laminates. The undamaged damping 
value was bv far the highest of the specimens tested (the [90Jg 
laminate was higher, but w*as only tested to obtain model 
constants), and was closely predicted by the mathematical 
model as seen in Figure 9. The plot of the logarithmic 
decrement vs. S* this laminate appears in Figure 35, with the 
corresponding model predictions. The initial and final 
predictions corres p ond very well to the measured values. The 
measured value at the middle damage state is slightly higher 
than the final measurement, and is attributed to experimental 
error. 
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Laminate H. 

Tne (O/P'DjOj (equivalently 1 0/90/0/ 90) s ) stacking sequence 
is similar to that of the I0/9D/0] S except for the addition of 


two 90 degree plies adjacent to the laminate mid-plane. Tne 
damping for the undamaged |0/90H« specimens was found to be 
almost identical in that measured for the 1 0/90/ 01* specimens, 
and was closely predicted by the model. The damage 
progressed vcry‘r*pid)y for this particular sucking sequence, 
in the 90 degree plies. This laminate began to show matrix 
cracking at 77 *> ksi and quickly reached a matrix crack 
saturation level at a stress level of 85 ksi, and failed at 875 
ksi. The damping increased dramatically in these specimens, 
and is seen in Figure 16. This increase is predicted very well by 
the model 
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Laminate r. f9Q/Q/9Cri h 

Tne sucking sequence for laminate F was similar to the 
I0/905) s series, except for one of the 90 degree plies being 
• outside the 0 degree ply. This difference has no effect on axial 
stiffness, but it does affect the bending stiffness. Tne damping 
of the undamaged [90/0/902)5 specimen is seen to be 
significantly higher than the [0/903)5 in Figure 9. For the 
190/0/902)5 sucking sequence, the matrix cracks gTew ven^ 
quickly in the outer 90 degree plies, and more slowly in the 
inner 90 degree plies. Tnis laminate showed a very large 
increase in damping with increased damage, which was 
closely predicted by the model. Tnis resulting function is 
plotted in Figure 37. 

Laminate Comparisons 

Tne results which have been presented indicate a strong 
dependence of damping on laminate stacking sequence and 
matrix cracking. It was postulated in the model development 
that the damping of a laminate would be greatly dependent 
-upon the 90 degree plies, and only slightly dependent upon the 
' v 0 degree plies. Tnis assumption is supported by the fact that 
the camping of the [90)g laminate was found to be 7 times 
greater than the damping of the [Ojg laminate. However, the 
extensive testing of cross-plv laminates indicated that the 
percentage of 90 degree plies in a specimen has very little 
effect on the material camping of the laminate. Figure IS is a 
plot of the measured logarithmic decrement for all of the 
undamaged specimens against the corresponding percentage of 
90 degree plies in the specimen. For example, the [0/90/0] s 
laminate contains 3333 percent 90 degree plies. No trend can 
be seen in this plot which would indicate that the damping 
changes with the percentage of 90 degree plies. 
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shown in Figure 19, ind indioitt interesting function of the 
k*c*rithnruc decrement with respect to the number and locator 
of 90 degree plies in • laminate. The trend appears to be 
exponential in lorm, although since no exact correlation can be* 
made to accoum lor this function, a curve hi was not applied. 
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The major factor contributing to the relative damping 
magnitudes for the undamaged specimens was seen to be the 
placement of the 90 degree plies with respect to the mid- 
plane. . A plot similar to Figure IS was created which 
weighted each 90 degree ply by the distance from the mid- 
piane, and related this value to the measured logarithmic 
decrement. The value used for the domain is a weighted 
percentage of 90 degree plies given bv the following formula: 

J 


~* 5 > 


where J is the number of 90 degree plies and N is the total 
number of piies in the laminate. A laminate with a significant 
fraction of the 90 degree plies far 2 way from the mid-plane 
will have a high value of 'F, while a laminate with the same 
percentage of 90 degree piies closer to the mid-plane, will 
have a low value of V F. The graph of this relationship is 


Fcramufte of 9C Tiim In Uuniniu 

c 

Fijrun 1$. MunxmC Lopxrilhcruc !D*cr«ar*ent v*. Weghiad * 91 The* 


Figure 20 is a plot of the increase in the logarithmic 
decrement for each stacking sequence, from the undamaged 
state, to the matrix crack saturation damage state. The 
corresponding increase predicted by the model is also shown 
for each specimen. In all of the specimens except the 1 903/0)5 
specimen, the predicted increase is lower than the measured 
increase. This difference is possibly due to the presence of 
additional damage which may have been introduced during 
the loading process. The stacking sequences which were 
greatly underpredicted, I90/0/90) s (23 percentage points 
under), and the 1 93/0/ 90)5 (25.7 perc en tage points under), both 
had 90 degree plies on the outer layer. As the number of 
. matrix cracks directly under the strain gage increased for both 
of these laminates, the movement of the crack faces against 
the strain gage resulted in an increased damping measurement. 
This did not occur in the [903/0)5 laminate because the 
damping was not measured at a damage level vvhich contained 
a great number of matrix cracks. From the damage-dependent 
damping plots of the individual laminates given throughout 
the results section, it cam be seen that for these specimens with 
the 90 degree plies at the outer surface, the model makes 
accurate predictions at earlier matrix crack damage states. 
For specimens with 90 degree plies at the outer surface 
containing a large number of matrix cracks, the model does not 
seem to correspond well with the experimental measurements. 
For all of the laminates with 0 degree plies at the outer 
surface, the model provides reasonable predictions of the 
increase in damping at all matrix crack levels. 


The material damping of several graphite/epoxy cross- 
ply laminates was measured as a function of matrix cracking. 
The testing configuration was a cantilevered beam, 
mechanically excited to vibrate in the first mode. The free 
response was captured, and the logarithmic decrement was 
calculated from the peak decay of the strain amplitude. Six 
different laminates were tested at several levels of matrix 
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crmck (Unuge. A mathem* oca! mode] was used to predict the 
damping at *11 levels of damage* including the undamaged 
s p ecimens. 

The damping values lor all of the undamaged specimens 
were predicted fry the mamemanca) model. The percent error 
lor this set of measurements and predictions was less than i 
13% for all of the specimens tested. The measured increases 
were very significant in all of the specimens. It was seen that 
matrix cracks in 90 degree plies which were far away from the 
mid -plane caused a greater increase in laminate damping 
values* than matrix cracks near the mid-plane. 
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8.4 A Model for Predicting Damage Dependent Damping in 
Laminated Composites 
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A MODEL FOR PREDICTING DAMAGE DEPENDENT 
DAMPING IN LAMINATED COMPOSITES 

S. Kalyanasundaram and D. H. Allen 

Aerospace Engineering Department, Texas AAiM University, College Station, TX 77843, U.S.A. 


Abstract — A constitutive model has been developed for predicting the dynamic properties of 
composite beams which accounts for the influence of microstructural damage. The model assumes the 
damage state to be fixed. The model can account for both the damping increase and stiffness loss due 
to microstructural damage. Damage dependent material constants obtained from a baseline data set 
can be used to predict damping changes for several composite layups. The validity- of the mode! has 
being verified by comparing the experimental results to model predictions for crossply laminates with 
a variety of stacking sequences. 


NOTATION 

A cross sectional area of the beam 

£ Young's modulus of elasticity 

5 total number of plies in the laminate 

s total surface are 2 of matrix cracks 

V maximum strain energy of deflected beam 

H* deflection of the beam 

6 logarithmic decrement 

P mass densiry/unit length 

l length cf the beam 

P complex curvature 

frequency parameter 
Tj loss factor 

compicx frequency 

£5 complex cnensional modulus of the central ply 
area moment of inertia of the cross section 
£* complex Young's modulus 

£' storage modulus 

£" loss modulus 

A f complex bending moment 

E~, complex exiensional modulus of /th plv 

£", extensional storage modulus of the jth ply 
£" extensional loss modulus of the ;th ply 
£L complex modulus of C c plies 

£t complex modulus of 90 c plies 

tan flexural loss tangent 

INTRODUCTION 

The process of ultimate failure of advanced com- 
posite materials is preceded by a sequence of 
microstructural and macrostructural events termed 
as damage. These events may be due to transverse 
matrix cracking, delamination, fiber breaking and 
fiber matrix debondingfl— 4]. This damage results 
in stiffness loss which in turn may have a significant 
effect on the dynamic response of structural 
components[5j. Experimental results indicate that 
material damping may be more sensitive to damage 
development than stiffness loss. Damping in the 
presence of damage has been observed as high as 
350% greater than the reference undamaged 
camping although the corresponding reduction in 
resonant frequencies is only about 5%. 

Schult2 and Warwick [6] used oscillatory flexural 


loads of zero mean value to fatigue cantilever 
beam specimens of scotchply laminates. During 
this fatiguing, the vibration response of the speci- 
mens was measured to determine if the response 
chances were indicative of fatigue. It was found 
that the increase in damping correlated fairly well 
with damage. For specimens which exhibited 
significant crack damage the damping ratio 
increased by 100%. In a closely related in- 
vestigation [7] similar results were reported for os- 
cillatory tensile loads of non-zero mean value. A 
torsional pendulum was used to measure the tor- 
sional storage modulus and damping capacity- of 
carbon fiber reinforced bars before and after 
cracks were introduced through static torsion 
testing [S]. The specific damping capacity' increased 
by nearly 80% after the introduction of damage. 
Chandra ti aL[ 9] measured damping in crossply 
giass/epoxv composite beams before and after the 
introduction of damage due to static and fatigue 
loads. It was found that there was a considerable 
increase in damping with the increase in number of 
cycles of fatigue loading and the level of prestress- 
ing. For the fatigue load the reduction in storage 
modulus was very small after five cycles, but the 
damping showed a monotonic increase with in- 
creasing cycles. The damping was found to be in- 
dependent of amplitude of oscillation up to about 
0.04% strain. 

Plunkett [10] measured damping factors of can- 
tilever beams made of crossply glass/epoxy 
laminates as they were vibrated at different strain 
levels in the firs: and second bending mode simul- 
taneously. It was found that the damping factor 
was substantially increased by large strain induced 
matrix cracking in laminae with transverse fibers. 
The damping factor was independent of current 
strain history or mode shape and was a function of 
damace. It was concluded that this behavior makes 
it possible to use a linear superposition method for 


13 3 


114 


S. Kalyanasundaxam and 0. R Aluen 


dynamic analysis. Smith et a/. [11] have measured 
maieriaJ damping of laminated graphite/epoxy 
specimens as a function of matrix cracking. It was 
found that damping increased significantly lor 
various stacking sequences due to transverse 
matrix cracks. Alien et aL[l 2] have developed a 
theoretical model for predicting the damage 
dependent damping based on the theory of internal 
state variables to account for damage. Weits- 
man[13] has developed a continuum damage 
mode) for viscoelastic material. The special cases 
of uniaxial damage under uniaxial stress and inter- 
action of damage with moisture diffusion have 
been considered. 

In this paper a mode] is developed for predicting 
the structural dynamic response of composite 
beams in the presence of transverse matrix cracks. 
The validity of the theoretical formulation is 
verified by comparing the theoretical results with 
the experimental results obtained in [13] for 
different stacking sequences. 


mathematical model 

The determination of damping properties of 
beams usually employs the Euler-Bemoulli 
assumption. In this theory’ it is assumed that the 
beam element is in pure bending and shear effects 
arc neglecied. For highly anisotropic materials and 
in the presence of interply delaminations, the 
dynamic properties obtained on the basis of 
the Euler-Bemoulli beam theory’ may exhibit 
significant errors. Since the current study will con- 
centrate on crossplv laminated beams with trans- 
verse matrix cracks, the Euler-Bemoulli beam 
theory is assumed to be adequate for obtaining the 
dynamic properties. The governing differential 
equation for the dynamic response of an elastic 
Euler-Bemoulli beam is given by 


C* K’ K* 

^aF +M a?- a 


0) 


The solution for free vibrations can be written in 
the form 


*’„(*, t) = exp(iotf). (2) 

For sufficiently small damping, free vibrations are 
also approximately harmonic and from the cor- 
respondence principle of viscoelasticity’, the elastic 
Young’s modulus may be replaced by the cor- 
responding complex modulus, E where 

Er-E'+iEr = E'{\ + h) (3) 

and the flexural damping is given by 


tan 6 e 


The complex frequency equation is given by 


A 4 .’ 


pAVu>\ 


( 5 ) 


where A„ is given by the boundary’ condition. Since 
E* is complex, the frequency values are also com- 
plex and can be written as 

*= + <a>I . (6) 

Separating the real and imaginary pans of eqn (5) 
results in 
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tan 4> e = 
tan = 




5 


"V 


( 7 ) 


(S) 


(9) 


where 5 is the logarithmic decrement obtained 
from a free vibration decay experiment. 

The uniaxial constitutive equation for matrix 
cracking of elastic-composites under isothermal 
conditions is modelled by[14, 15]. 

+ 00 ) 

where cr, is an internal state variable representing 
matrix cracking in each ply. The changes in dam- 
ping and stiffness due to mi crostru crural damage 
are history dependent and spatially variable. For a 
fixed damage state the dynamic properties are 
slowly varying functions of time so that the com- 
posite material behaves like a linear viscoelastic 
body. Under sinusoidal loading, experimental 
results[10] indicate that the principle of super- 
position is valid for a damage state that is fixed and 
subjected to small amplitude vibrations. Based on 
these observations it is proposed that eqn (10) can 
be modified for sinusoidal loading at fixed damage 
state to 


< 7 „ = E 9 i„ + ra lt 


(ID 


where 




(4) 


r = r + ir 
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Further insight into the effects of fixed damage 
on damping can be gained by studying the energy 
dissipated in a cycle of sinusoidal motion. There- 
fore, 


iW *= J 



(12) 

» d(Dt) 

(13) 


evidence that damping increases with micros truc- 
tural damage arid the resonant frequencies 
decrease. Equation (18) can therefore be written as 


ET + AE" 

’’"f+af' 


( 20 ) 


Application of the model to crossply laminates 

The moment-curvature relationship for a visco- 
elastic Euler-Bemoulli beam vibrating sinusoidal- 
ly with frequency u can be written as 


M = 1 
ETI p 


(21) 


In order to obtain a meaningful expression for 
energy dissipation the real part of and 
should be used in the evaluation of (13). Thus, 

rcal [S)] = "^ sin<1 " (14) 

Real(*„) = Real [(F- ifjc* e~ 

-f (T — //") a. A e^'] 

= [E* sin wt — F cos wr] e A 
-r[/' sin cut — r cos a)f]a, A . 

(15) 

Using (15) and (14) in (13) the energy dissipated 
per cycle is given by 

Aw = ‘rE* + 7r/ r a 1A e A . (16) 

Note that for a 1A = 0 (no damage) the above 
expression reduces to the case of energy dissipation 
in a linear viscoelastic body with no damage. 

The strain energy stored during the peak dis- 
placement is 

L/ = [F£-+J'^ a , A ]/2. (17) 

The loss factor is 

<1! > 

E e* + I 

Furthermore, it will be assumed that 

ra u = AF£ A (19a) 

Ta 1A =AFc a , (19b) 

where AF is the increase in loss modulus and A£ r 
is the decrease in storage modulus due to damage. 
This assumption is made based on the experimental 


For plies of uniform thickness the complex flexural 
modulus is given in terms of the complex ply 
moduli and stacking geometry by[16]: 

g r*n “I 

£* - -^5 [ I E l ~ 3/ + 1) j for even N 

■wtf + 

for odd N. 

In this work symmetric laminated beams composed 
solely of longitudinal and transverse plies will be 
considered. For this case eqn (22) reduces to 


E-= — [aE- L +bE’l 


(23) 


where a and b are functions of stacking geometry. 
The flexural damping is given by 


tan d>£ ~ Tj — 


aEl + bE" r 
aE' L +bE' T ' 


(24) 


In the presence of transverse matrix cracks in the 
90 c plies, it is assumed here that j E' r and E m r will be 
affected by damage and the values of E' L and El 
remain unchanged since no cracking occurs in the 
0° plies. 

' .Comparing eqns'(20) and (24), the loss factor in 
the presence of transverse matrix cracks is thus 
given by 


aE m L + CAEr+ *>Et 
V ~ aE' L + CAEV + bE'r 


(25) 


where CAEV and CAEV are functions of the total 
surface area of matrix cracks and the location of 
‘the cracks. AEV and AEV as functions of crack 
damage can be extracted from a single simplified 
stacking sequence. This information may then be 
utilized in (25) to predict the damping as a function 
of damage for different stacking sequences. 
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DISCUSSION OF RESULTS 

Experimental results [11] obtained from crossply 
laminates with different stacking sequences will be 
used to verify the validity of the model developed 
in the previous section. By testing laminates built 
from 0 or 90* plies, the values of E' L and E' T can 
be obtained from resonant frequency measure- 
ments through eqn (8). Using the damping 
measurements, the values of E m L and EV can be 
obtained from eqn (4). These values can be used to 
predict the damping of any crossply laminates with 
various combinations of 0 and 90* plies. Figure 1 is 
a comparison of experimental and theoretical 
damping values for different laminates. For 
polymer matrix composites the viscoelasticity of 
the matrix is the major source of energy dis- 
sipation. Hence, the material damping of a crossply 
laminated beam is dominated by the location and 
the number of 90° plies. The predicted and 
experimental damping values indicate that the 
flexural damping increases with distance from the 
midplane of the beam. 

For predicting the damage dependent damping 
for different stacking sequences, the increase in 
loss modulus as a function of matrix crack damage 
can be constructed by experimentally measuring 
the damping increase in a [0/90/0], laminate. The 
change in loss factor as a function of transverse 
matrix crack surface area is given experimentally 
by 

A Ty (s) « 0.0002(s) - 1.7 1 OS x 1 O' 5 s\ (26) 

For this laminate the damping increased by 35%, 
while the reduction in the resonant frequency was 
about 0.07%. Using eqn (25) for a [0/90/0], 
laminate, the change in the loss modulus as a func- 



rig. 1. Damping values for various laminates. 


tion of surface area of matrix cracks is given by 

A rAs) «= 57.97 x 10“ 4 [0.0002 (s) 

- 1.7108 x 10-*VJ. (27) 

The experimental change in the storage modulus is 
negligible. 

The experimental and theoretical predictions of 
damping of a [90/0/90], laminate at two different 
damage stales are illustrated in Fig. 2. The de- 
scriptions of the damage states are given in Table 
1. For this laminate there are more transverse 
matrix cracks at the outer 90* plies than the 90* 
plies near the midplane. The damping increases by 
nearly 100%. The agreement between theoretical 
and experimental results is good. The reduction in 
the first mode resonant frequency (an indication of 
stiffness loss) is only about 7%. This trend clearly 
indicates that damping changes are more sensitive 
to the microstructuraliv induced damping changes 
than the stiffness loss. Figure 3 depicts the effect of 
the matrix damage on a [0/90/0/90], laminate. The 
corresponding damage states arc given in Table 2. 
The theoretical model predicts the increase in 
damping due to transverse matrix cracks to be a 
function of both total matrix crack area and the 
location of the 90° ply in the laminate. It is interes- 
ting to note that for approximately the same 
amount of total matrix crack area the damping 
change in a [90/0/90], laminate is several times 
higher than the change ’in the [0/90/0/90], 
laminate. This is due to the fact that the matrix 
crack damage in 90* plies farther away from the 
midplane of the beam will lead to a larger increase 
in damping than the 90° plies near the midp) 2 ne. 
This trend is captured by both the theoretical and 
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Fig. 2. Damage dependent damping values for [90/0/90], 
laminate. 
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Fic. 3. Damage dependent damping values for [0/90 ] 2j 
laminate. 
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Fig. 4. Damage dependent damping values for 
[90/0/90;], laminate. 

experimental results. Figure 4 illustrates the effect 
of matrix cracks on a [90/0/90;], laminate (see 
Table 3). The agreement between the theoretical 
and experimental damping values for the damage 
states is good. However, the theoretical model 
tends to overpredict the increase in damping from 
undamaged damping values as compared to the 
experimental results. 

CONCLUSIONS 

A constitutive model has been developed for 
including the effect of microstructural damage on 
the dynamic response of laminated composites. 


Tabic 1. Description of damage state for 
(90/0/90), alminate 


Sucking 

geometry 

Damage state I 
(cracks/in.) 

Damage state II 
(cracks/in.) 

90 

7.9 

21.3 

0 

- 

- 

90 

0 

2.9 

90 

0 

2.9 

0 

- 

- 

90 

7.9 

30.2 


Table 2. Description of damage state for 
[0/90];, laminate 


Sucking 

geometry 

Damage sute I 
(cracks/in.) 

Damage sute 11 
(cracks/in.) 

0 

_ 

_ 

90 

5.0 

10.3 

0 

- 

- 

90 

6.3 

13.0 

90 

6.3 

13.0 

0 

- 

- 

90 

3.7 

6.7 

0 

- 

- 


Table 3. Description of damage state for 
[90/0/90J, laminate 


Stacking 

geometry’ 

Damage state I 
(cracks/in.) 

Damage state II 
(cracks/in.) 

90 

10.9 

24.7 

0 

- 

- 

90 

.3.9 

27.2 

90 

3.9 

27.2 

90 

3.9 


90 

3.9 

27.2 

0 

- 

- 

90 

11.9 

26.7 


This model is valid for a fixed damage state and 
small amplitude vibrations. Under these assump- 
tions the composite material behaves like a quasi- 
linear viscoelastic body. The validity of the 
theoretical formulation has been verified for 
crossply graphite-epoxv laminates with a variety' of 
stacking sequences. Damage dependent material 
constants obtained from [0/90/0], laminate have 
been used to predict the increase in damping in 
[90/0/90],, [0/90/0/90], and [90/0/90;], laminates 
for different damage states. Damping is found to 
be more sensitive to microstructural damage than 
the stiffness loss. Thus, damping holds some prom- 
ise for further studying the damage development in 
composite structural components. 
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